basic input variables¶
This document lists and provides the description of the name (keywords) of the basic input variables to be used in the input file for the abinit executable.
accuracy¶
Mnemonics: ACCURACY
Mentioned in topic(s): topic_Planewaves, topic_SCFControl
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [17/1121] in all abinit tests, [3/156] in abinit tutorials
Allows to tune the accuracy of a groundstate or DFPT calculation optdriver=0 or 1, by setting automatically the variables according to the following table:
accuracy  1  2  3  4  5  6 

ecut  E_min  E_med  E_med  E_max  E_max  E_max 
pawecutdg  ecut  ecut  1.2 * ecut  1.5 * ecut  2 * ecut  2 * ecut 
fband  0.5  0.5  0.5  0.5  0.75  0.75 
boxcutmin  1.5  1.8  1.8  2.0  2.0  2.0 
bxctmindg  1.5  1.8  1.8  2.0  2.0  2.0 
pawxcdev  1  1  1  1  2  2 
pawmixdg  0  0  0  0  1  1 
pawovlp  10  7  7  5  5  5 
pawnhatxc  0  1  1  1  1  1 
tolvrs  1.0d3  1.0d5  1.0d7  1.0d9  1.0d10  1.0d12 
tolmxf  1.0d3  5.0d4  1.0d4  5.0d5  1.0d6  1.0d6 
optforces  1  1  2  2  2  2 
timopt  0  0  1  1  1  1 
npulayit  4  7  7  7  15  15 
nstep  30  30  30  30  50  50 
prteig  0  0  1  1  1  1 
prtden  0  0  1  1  1  1 
accuracy = 4 corresponds to the default tuning of ABINIT. It is already a very accurate tuning. For a parallel calculation, timopt is enforced to be 0. E_min, E_med and E_max may be read from the pseudopotential file (available only for XML PAW atomic data files). If E_min, E_med and E_max are not given in the pseudopotential file, ecut must be given in the input file and E_max=E_med=E_max=ecut. If the user wants to modify one of the input variable automatically tuned by accuracy, they must put it in the input file. The other input variables automatically tuned by accuracy will not be affected. accuracy = 0 means that this input variable is deactivated.
For the other values of optdriver, many of the above input variables have no meaning, so the accuracy has to be tuned by the user (e.g. for GW calculations, perform convergence studies with respect to ecuteps and other relevant input variables).
acell¶
Mnemonics: CELL lattice vector scaling
Characteristics: EVOLVING, LENGTH
Mentioned in topic(s): topic_UnitCell
Variable type: real
Dimensions: (3)
Commentdims: represented internally as acell(3,nimage)
Default value: 3 * 1
Added in version: before_v9
Test list (click to open). Very frequently used, [1094/1121] in all abinit tests, [144/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives the length scales by which dimensionless primitive translations (rprim) are to be multiplied. By default, given in Bohr atomic units (1 Bohr=0.5291772108 Angstroms), although Angstrom can be specified, if preferred, since acell has the LENGTH characteristics. See further description of acell related to the rprim input variable, the scalecart input variable, and the associated internal %rprimd input variable.
Note that acell is NOT the length of the conventional orthogonal basis vectors, but the scaling factors of the primitive vectors. Use scalecart to scale the cartesian coordinates.
angdeg¶
Mnemonics: ANGles in DEGrees
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_UnitCell
Variable type: real
Dimensions: (3)
Default value: None
Comment: deduced from ‘rprim‘
Added in version: before_v9
Test list (click to open). Moderately used, [61/1121] in all abinit tests, [0/156] in abinit tutorials
Gives the angles between directions of primitive vectors of the unit cell (in degrees), as an alternative to the input array rprim. Will be used to set up rprim, that, together with the array acell, will be used to define the primitive vectors.
 angdeg(1) is the angle between the 2^{nd} and 3^{rd} vectors,
 angdeg(2) is the angle between the 1^{st} and 3^{rd} vectors,
 angdeg(3) is the angle between the 1^{st} and 2^{nd} vectors,
If the three angles are equal within 1.0d12 (except if they are exactly 90 degrees), the three primitive vectors are chosen so that the trigonal symmetry that exchange them is along the z cartesian axis:
R1 = ( a, 0, c)
R2 = (a/2, sqrt(3)/2*a, c)
R3 = (a/2, sqrt(3)/2*a, c)
where a 2 + c 2 = 1.0 If the angles are not all equal (or if they are all 90 degrees), one will have the following generic form:
 R1 = (1, 0, 0)
 R2 = (a, b, 0)
 R3 = (c, d, e)
where each of the vectors is normalized, and form the desired angles with the others.
ecut¶
Mnemonics: Energy CUToff
Characteristics: ENERGY
Mentioned in topic(s): topic_Planewaves
Variable type: real
Dimensions: scalar
Default value: None
Added in version: before_v9
Test list (click to open). Very frequently used, [1121/1121] in all abinit tests, [156/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Used to define the kinetic energy cutoff which controls the number of planewaves at given k point. The allowed plane waves are those with kinetic energy lower than ecut, which translates to the following constraint on the planewave vector \vec{G} in reciprocal space
All planewaves inside this basis sphere centered at k are included in the basis (except if dilatmx is defined). The cutoff can be specified in Ha units (the default), Ry, eV or Kelvin, since ecut has the ENERGY characteristics. (1 Ha = 27.2113845 eV)
This is the single parameter which can have an enormous effect on the quality of a calculation; basically the larger ecut is, the better converged the calculation is. For fixed geometry, the total energy MUST always decrease as ecut is raised because of the variational nature of the problem.
Usually one runs at least several calculations at various ecut to investigate the convergence needed for reliable results.
For kpoints whose coordinates are build from 0 or ½, the implementation of timereversal symmetry that links coefficients of the wavefunctions in reciprocal space has been realized. See the input variable istwfk. If activated (which corresponds to the default mode), this input variable istwfk will allow to divide the number of plane wave (npw) treated explicitly by a factor of two. Still, the final result should be identical with the ‘full’ set of plane waves.
See the input variable ecutsm, for the smoothing of the kinetic energy, needed to optimize unit cell parameters.
einterp¶
Mnemonics: Electron bands INTERPolation
Mentioned in topic(s): topic_ElecBandStructure, topic_SelfEnergy
Variable type: real
Dimensions: (4)
Default value: [0, 0, 0, 0]
Added in version: before_v9
Test list (click to open). Rarely used, [5/1121] in all abinit tests, [1/156] in abinit tutorials
 libxc: t42.abi
 tutorespfn: teph4isotc_3.abi
 v8: t04.abi
 v9: t57.abi, t62.abi
This variable activates the interpolation of the electronic eigenvalues. It can be used to interpolate KS eigenvalues at the end of the GS run or to interpolate GW energies in sigma calculations (optdriver = 4). The kpath can be specified with kptbounds and nkpath. einterp consists of 4 entries. The first element specifies the interpolation method.
 0 → No interpolation (default)
 1 → Starfunction interpolation (modified ShanklandKoellingWood Fourier interpolation scheme, see [Pickett1988]).
The meaning of the other entries depend on the interpolation technique selected. In the case of starfunction interpolation:
 einterp(2): Number of starfunctions per abinitio kpoint
 einterp(3): If nonzero, activate Fourier filtering according to Eq. 9 of [Uehara2000]. In this case, rcut is given by einterp(2) * Rmax where Rmax is the maximum length of the lattice vectors included in the star expansion
 einterp(4): Used if einterp(3) /= 0. It defines rsigma in Eq. 9
iscf¶
Mnemonics: Integer for SelfConsistentField cycles
Mentioned in topic(s): topic_SCFAlgorithms, topic_TDDFT, topic_ElecBandStructure
Variable type: integer
Dimensions: scalar
Default value: 17 if %usepaw == 1,
0 if usewvl == 1,
7 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [335/1121] in all abinit tests, [51/156] in abinit tutorials
 bigdft: t00.abi, t02.abi, t03.abi …
 bigdft_paral: t01.abi, t01.abi …
 builtin: testin_bigdft.abi …
 etsf_io: t02.abi, t30.abi …
 fast: t09.abi, t11.abi, t12.abi …
 gpu: t01.abi …
 libxc: t13.abi, t19.abi, t81.abi …
 mpiio: t22.abi, t62.abi, t62.abi …
 paral: t05.abi, t05.abi, t05.abi …
 psml: t03.abi, t04.abi, t07.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdmft_1.abi, tgswvl_1.abi, tgswvl_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase3_5.abi, tbs_1.abi, tgw1_1.abi …
 v1: t44.abi, t45.abi, t46.abi …
 v2: t04.abi, t06.abi, t07.abi …
 v3: t09.abi, t13.abi, t16.abi …
 v4: t02.abi, t07.abi, t20.abi …
 v5: t02.abi, t05.abi, t08.abi …
 v6: t11.abi, t14.abi, t20.abi …
 v67mbpt: t04.abi, t11.abi, t14.abi …
 v7: t07.abi, t41.abi, t43.abi …
 v8: t04.abi, t07.abi, t12.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t03.abi, t11.abi, t12.abi …
Controls the selfconsistency algorithm.
Positive values correspond to the usual choice for doing the usual ground state (GS) calculations or for structural relaxations, where the potential has to be determined selfconsistently while negative values correspond to nonselfconsistent calculations.
Different SCF algorithms (> 0) can be selected:
 0 → SCF cycle, direct minimization scheme on the gradient of the wavefunctions. This algorithm is faster than diagonalisation and mixing but is working only for systems with a gap. It is implemented only on the wavelet basis set, when usewvl = 1.

1 → get the largest eigenvalue of the SCF cycle (DEVELOP option, used with irdwfk = 1 or irdwfq = 1)

2 → SCF cycle, simple mixing of the potential with the preconditioned potential residual (in the usual case, the latter is defined by diemix, diemac and dielng)
 3 → SCF cycle, Anderson mixing of the potential
 4 → SCF cycle, Anderson mixing of the potential based on the two previous iterations
 5 → SCF cycle, CG based on the minim. of the energy with respect to the potential
 6 → SCF cycle, CG based on the minim. of the energy with respect to the potential (alternate algo., DEVELOP)
 7 → SCF cycle, Pulay mixing of the potential based on the npulayit previous iterations
 12 → SCF cycle, simple mixing of the density
 13 → SCF cycle, Anderson mixing of the density
 14 → SCF cycle, Anderson mixing of the density based on the two previous iterations
 15 → SCF cycle, CG based on the minim. of the energy with respect to the density
 16 → SCF cycle, CG based on the minim. of the energy with respect to the potential (alternate algo., DEVELOP)
 17 → SCF cycle, Pulay mixing of the density based on the npulayit previous iterations
Warning
Other positive values, including zero ones, are not allowed.
Such algorithms for treating the “SCF iteration history” should be coupled with accompanying algorithms for the SCF “preconditioning”. See the input variable iprcel. The default value iprcel = 0 is often a good choice, but for inhomogeneous systems, you might gain a lot with iprcel = 45.
(Warning: if iscf > 10, at present (v4.6), the energy printed at each SCF cycle is not variational  this should not affect the other properties, and at convergence, all values are OK)
 In the normconserving case, the default option is iscf = 7, which is a compromise between speed and reliability. The value iscf = 2 is safer but slower.
 In the PAW case, default option is iscf = 17. In PAW you have the possibility to mix density/potential on the fine or coarse FFT grid (see pawmixdg).
 Note that a Pulay mixing (iscf = 7 or 17) with npulayit = 1 (resp. 2) is equivalent to an Anderson mixing with iscf = 3 or 13 (resp. 4 or 14).
 Also note that:
 when mixing is done on potential (iscf < 10), total energy is computed by “direct” decomposition.
 when mixing is done on density (iscf >= 10), total energy is computed by “double counting” decomposition. “Direct” and “double counting” decomposition of energy are equal when SCF cycle is converged. Note that, when using GGA XC functionals, these decompositions of energy can be slightly different due to imprecise computation of density gradients on FFT grid (difference decreases as size of FFT grid increases  see ecut for NC pseudopotentials, pawecutdg for PAW).
Other (negative) options:

2 → a nonselfconsistent calculation is to be done; in this case an electron density rho(r) on a real space grid (produced in a previous calculation) will be read from a disk file (automatically if ndtset = 0, or according to the value of getden if ndtset/=0). The name of th density file must be given as indicated here. iscf = 2 would be used for band structure calculations, to permit computation of the eigenvalues of occupied and unoccupied states at arbitrary k points in the fixed self consistent potential produced by some integration grid of k points. Due to this typical use, ABINIT insist that either prtvol > 2 or prteig does not vanish when there are more than 50 k points. To compute the eigenvalues (and wavefunctions) of unoccupied states in a separate (nonselfconsistent) run, the user should save the selfconsistent rho(r) and then run iscf = 2 for the intended set of kpoints and bands. To prepare a run with iscf = 2, a density file can be produced using the parameter prtden (see its description). When a selfconsistent set of wavefunctions is already available, abinit can be used with nstep = 0 and the adequate value of prtden, see tests/v2/Input/t47.abi .

3 → like 2, but initialize occ and wtk, directly or indirectly (using ngkpt or kptrlatt) depending on the value of occopt. For GS, this option might be used to generate Densityofstates (thanks to prtdos), or to produce STM charge density map (thanks to prtstm). For RF, this option is needed to compute the response to ddk perturbation.

1 → like 2, but the nonselfconsistent calculation is followed by the determination of excited states within TDDFT. This is only possible for nkpt = 1, with kpt = 0 0 0. Note that the oscillator strength needs to be defined with respect to an origin of coordinate, thanks to the input variable boxcenter. The maximal number of KohnSham excitations to be used to build the excited state TDDFT matrix can be defined by td_mexcit, or indirectly by the maximum KohnSham excitation energy td_maxene. Only “LDAlike” XC kernels are implemented, from native ABINIT XC functionals with specific values of ixc=[0,1,7,8,20,21,22]. In case you use a XC functional from the libxc, corresponding to these values, there is no automatic translation, so that you have to specify explicitly the ixc value among the above ones. As an example, if the ixc from your pseudopotentiel is 1012, this corresponds to the same functional (LDA  Perdew Wang 1992) as ixc=7, and you should specify ixc=7.
ixc¶
Mnemonics: Index of eXchangeCorrelation functional
Mentioned in topic(s): topic_xc, topic_Hybrids, topic_TDDFT
Variable type: integer
Dimensions: scalar
Default value: 1
Comment: Default corresponds to Teter parametrization. However, if all the pseudopotentials have the same value of pspxc, the initial value of ixc will be that common value
Added in version: before_v9
Test list (click to open). Moderately used, [389/1121] in all abinit tests, [24/156] in abinit tutorials
 bigdft: t00.abi, t06.abi, t07.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_bigdft.abi, testin_etsf_io.abi, testin_libxc.abi …
 etsf_io: t00.abi, t09.abi …
 gpu: t01.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_82.abi, tsv3_03.abi, tsv3_04.abi …
 tutorespfn: teph4mob_1.abi, teph4mob_4.abi, tlw_1.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t10.abi, t11.abi, t12.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t06.abi, t07.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t07.abi, t21.abi, t23.abi …
 v6: t01.abi, t02.abi, t12.abi …
 v67mbpt: t09.abi, t50.abi, t51.abi …
 v7: t11.abi, t22.abi, t26.abi …
 v8: t21.abi, t22.abi, t32.abi …
 v9: t17.abi, t19.abi, t31.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Controls the choice of exchange and correlation (xc). The list of XC functionals is given below. Positive values are for ABINIT native library of XC functionals, while negative values are for calling the much wider set of functionals from the ETSF LibXC library (by M. Marques), available at the LibXC home page Note that the choice made here should preferably agree with the choice made in generating the original pseudopotential, except for ixc = 0 (usually only used for debugging). A warning is issued if this is not the case. Unfortunately, pseudopotential (or PAW) generators for hybrid functionals and mGGA are currently under development, so that one usually uses GGA or LDA pseudopotentials instead. The error should be limited when GGA or LDA pseudopotentials with semicore states are used. Still this is a noncontrolled error. Moreover, the choices ixc = 1, 2, 3 and 7 are fits to the same data, from CeperleyAlder, and are rather similar, at least for spinunpolarized systems. The choice between the nonspinpolarized and spinpolarized case is governed by the value of nsppol (see below).
Native ABINIT XC functionals
NOTE: in the implementation of the spindependence of these functionals, and in order to avoid divergences in their derivatives, the interpolating function between spinunpolarized and fullyspinpolarized function has been slightly modified, by including a zeta rescaled by 1.d01.d6. This should affect total energy at the level of 1.d6Ha, and should have an even smaller effect on differences of energies, or derivatives. The value ixc = 10 is used internally: gives the difference between ixc = 7 and ixc = 9, for use with an accurate RPA correlation energy.

0 → NO xc.

1 → LDA or LSD, Teter Pade parametrization (4/93, published in [Goedecker1996], which reproduces PerdewWang 92 [Perdew1992a] (which reproduces CeperleyAlder [Ceperley1980]!).
 2 → LDA, PerdewZungerCeperleyAlder (no spinpolarization) [Perdew1981]
 3 → LDA, old Teter rational polynomial parametrization (4/91) fit to CeperleyAlder data (no spinpolarization) [Ceperley1980]
 4 → LDA, Wigner functional (no spinpolarization)
 5 → LDA, HedinLundqvist functional (no spinpolarization) [Hedin1971]
 6 → LDA, “Xalpha” functional (no spinpolarization)
 7 → LDA or LSD, PerdewWang 92 functional [Perdew1992a]
 8 → LDA or LSD, xonly part of the PerdewWang 92 functional [Perdew1992a]

9 → LDA or LSD, x and RPA correlation part of the PerdewWang 92 functional [Perdew1992a]

11 → GGA, PerdewBurkeErnzerhof GGA functional [Perdew1996]
 12 → GGA, xonly part of PerdewBurkeErnzerhof GGA functional [Perdew1996]
 13 → GGA potential of van LeeuwenBaerends [VanLeeuwen1994], while for energy, PerdewWang 92 functional [Perdew1992a]
 14 → GGA, revPBE of [Zhang1998]
 15 → GGA, RPBE of [Hammer1999]
 16 → GGA, HTCH93 of [Hamprecht1998]
 17 → GGA, HTCH120 of [Boese2000]  The usual HCTH functional.
 18 → (NOT AVAILABLE: used internally for GGA BLYP pseudopotentials from [Krack2005], available from the CP2K repository  use the LibXC instead, with ixc = 106131.

19 → (NOT AVAILABLE: used internally for GGA BP86 pseudopotentials from [Krack2005], available from the CP2K repository  use the LibXC instead, with ixc = 106132.

20 → FermiAmaldi xc ( 1/N Hartree energy, where N is the number of electrons per cell; G=0 is not taken into account however), for TDDFT tests. No spinpol. Does not work for RF.
 21 → same as 20, except that the xckernel is the LDA ( ixc = 1) one, for TDDFT tests.
 22 → same as 20, except that the xckernel is the BurkePetersilkaGross hybrid, for TDDFT tests.
 23 → GGA of [Wu2006].
 24 → GGA, C09x exchange of [Cooper2010].
 26 → GGA, HTCH147 of [Boese2000].
 27 → GGA, HTCH407 of [Boese2001].

28 → (NOT AVAILABLE: used internally for GGA OLYP pseudopotentials from [Krack2005], available from the CP2K repository  use the LibXC instead, with ixc = 110131.

40 → HartreeFock
 41 → PBE0, [Adamo1999].
 42 → PBE0⅓, [Guido2013].
ETSF Lib XC functionals
Note that you must compile ABINIT with the LibXC plugin in order to be able to access these functionals. The LibXC functionals are accessed by negative values of ixc . The LibXC contains functional forms for either exchangeonly functionals, correlationonly functionals, or combined exchange and correlation functionals. Each of them is to be specified by a threedigit number. In case of a combined exchange and correlation functional, only one such threedigit number has to be specified as value of ixc , with a minus sign (to indicate that it comes from the LibXC). In the case of separate exchange functional (let us represent its identifier by XXX) and correlation functional (let us represent its identified by CCC), a sixdigit number will have to be specified for ixc , by concatenation, be it XXXCCC or CCCXXX. As an example, ixc = 1012 gives the PerdewWang 1992 LDA ( ixc =7 as well), ixc = 020 gives the Teter93 LDA ( ixc =1 as well), while ixc = 101130 gives the PBE GGA ( ixc =11 as well). Note that for a metaGGA, the kinetic energy density is needed. This means having usekden = 1.
(S)LDA functionals (do not forget to add a minus sign, as discussed above)
 001 → XC_LDA_X [Dirac1930], [Bloch1929]
 002 → XC_LDA_C_WIGNER Wigner parametrization [Wigner1938]
 003 → XC_LDA_C_RPA Random Phase Approximation [GellMann1957]
 004 → XC_LDA_C_HL Hedin & Lundqvist [Hedin1971]
 005 → XC_LDA_C_GL Gunnarson & Lundqvist [Gunnarsson1976]
 006 → XC_LDA_C_XALPHA Slaters Xalpha
 007 → XC_LDA_C_VWN [Vosko1980]
 008 → XC_LDA_C_VWN_RPA [Vosko1980]
 009 → XC_LDA_C_PZ Perdew & Zunger [Perdew1981]
 010 → XC_LDA_C_PZ_MOD Perdew & Zunger (Modified) [Perdew1981] Modified to improve the matching between the low and high rs part
 011 → XC_LDA_C_OB_PZ Ortiz & Ballone (PZ) [Ortiz1994] [Ortiz1997] [Perdew1981]
 012 → XC_LDA_C_PW Perdew & Wang [Perdew1992a]
 013 → XC_LDA_C_PW_MOD Perdew & Wang (Modified) [Perdew1992a]; Added extra digits to some constants as in the PBE routine.
 014 → XC_LDA_C_OB_PW Ortiz & Ballone (PW) [Ortiz1994] [Ortiz1997] [Perdew1992a]
 017 → XC_LDA_C_vBH von Barth & Hedin [Barth1972]
 020 → XC_LDA_XC_TETER93 Teter 93 parametrization [Goedecker1996]
 022 → XC_LDA_C_ML1 Modified LSD (version 1) of Proynov and Salahub [Proynov1994]
 023 → XC_LDA_C_ML2 Modified LSD (version 2) of Proynov and Salahub [Proynov1994]
 024 → XC_LDA_C_GOMBAS Gombas parametrization [Gombas1967]
 025 → XC_LDA_C_PW_RPA Perdew & Wang fit of the RPA [Perdew1992a]
 027 → XC_LDA_C_RC04 RagotCortona [Ragot2004]
 028 → XC_LDA_C_VWN_1 Vosko, Wilk, & Nussair (1) [Vosko1980]
 029 → XC_LDA_C_VWN_2 Vosko, Wilk, & Nussair (2) [Vosko1980]
 030 → XC_LDA_C_VWN_3 Vosko, Wilk, & Nussair (3) [Vosko1980]
 031 → XC_LDA_C_VWN_4 Vosko, Wilk, & Nussair (4) [Vosko1980]
GGA functionals (do not forget to add a minus sign, as discussed above)
 084 → XC_GGA_C_OP_XALPHA oneparameter progressive functional (G96 version) [Tsuneda1999]
 085 → XC_GGA_C_OP_G96 oneparameter progressive functional (G96 version) [Tsuneda1999]
 086 → XC_GGA_C_OP_PBE oneparameter progressive functional (PBE version) [Tsuneda1999]
 087 → XC_GGA_C_OP_B88 oneparameter progressive functional (B88 version) [Tsuneda1999]
 088 → XC_GGA_C_FT97 Filatov & Thiel correlation [Filatov1997a] [Filatov1997]
Warning
this functional is not tested. Use at your own risks.
 089 → XC_GGA_C_SPBE PBE correlation to be used with the SSB exchange [Swart2009]
 090 → XC_GGA_X_SSB_SW Swart, Sola and Bickelhaupt correction to PBE [Swart2009a]
 091 → XC_GGA_X_SSB [Swart2009]
Warning
This functional gives NaN on IBM (XG20130608).
 092 → XC_GGA_X_SSB_D [Swart2009]
Warning
This functional gives NaN on IBM (XG20130608).
 093 → XC_GGA_XC_HCTH_407P HCTH/407+ [Boese2003]
 094 → XC_GGA_XC_HCTH_P76 HCTH p=7/6 [Menconi2001]
 095 → XC_GGA_XC_HCTH_P14 HCTH p=¼ [Menconi2001]
 096 → XC_GGA_XC_B97_GGA1 Becke 97 GGA1 [Cohen2000]
 097 → XC_GGA_XC_HCTH_A HCTHA [Hamprecht1998]
 098 → XC_GGA_X_BPCCAC BPCCAC (GRAC for the energy) [Bremond2012]
 099 → XC_GGA_C_REVTCA Tognetti, Cortona, Adamo (revised) [Tognetti2008]
 100 → XC_GGA_C_TCA Tognetti, Cortona, Adamo [Tognetti2008a]
 101 → XC_GGA_X_PBE Perdew, Burke & Ernzerhof exchange [Perdew1996] [Perdew1997]
 102 → XC_GGA_X_PBE_R Perdew, Burke & Ernzerhof exchange (revised) [Zhang1998]
 103 → XC_GGA_X_B86 Becke 86 Xalfa,beta,gamma [Becke1986]
 104 → XC_GGA_X_HERMAN Herman Xalphabeta GGA [Herman1969] [Herman2009]
 105 → XC_GGA_X_B86_MGC Becke 86 Xalfa,beta,gamma (with mod. grad. correction) [Becke1986] [Becke1986a]
 106 → XC_GGA_X_B88 Becke 88 [Becke1988]
 107 → XC_GGA_X_G96 Gill 96 [Gill1996]
 108 → XC_GGA_X_PW86 Perdew & Wang 86 [Perdew1986a]
 109 → XC_GGA_X_PW91 Perdew & Wang 91 [JP Perdew, in Proceedings of the 21^{st} Annual International Symposium on the Electronic Structure of Solids, ed. by P Ziesche and H Eschrig (Akademie Verlag, Berlin, 1991), p. 11. ] [Perdew1992] [Perdew1993]
 110 → XC_GGA_X_OPTX Handy & Cohen OPTX 01 [Handy2001]
 111 → XC_GGA_X_DK87_R1 dePristo & Kress 87 (version R1) [DePristo1987]
 112 → XC_GGA_X_DK87_R2 dePristo & Kress 87 (version R2) [DePristo1987]
 113 → XC_GGA_X_LG93 Lacks & Gordon 93 [Lacks1993]
 114 → XC_GGA_X_FT97_A Filatov & Thiel 97 (version A) [Filatov1997a]
 115 → XC_GGA_X_FT97_B Filatov & Thiel 97 (version B) [Filatov1997a]
 116 → XC_GGA_X_PBE_SOL Perdew, Burke & Ernzerhof exchange (solids) [Perdew2008]
 117 → XC_GGA_X_RPBE Hammer, Hansen & Norskov (PBElike) [Hammer1999]
 118 → XC_GGA_X_WC Wu & Cohen [Wu2006]
 119 → XC_GGA_X_mPW91 Modified form of PW91 by Adamo & Barone [Adamo1998]
 120 → XC_GGA_X_AM05 Armiento & Mattsson 05 exchange [Armiento2005] [Mattsson2008]
 121 → XC_GGA_X_PBEA Madsen (PBElike) [Madsen2007]
 122 → XC_GGA_X_MPBE Adamo & Barone modification to PBE [Adamo2002]
 123 → XC_GGA_X_XPBE xPBE reparametrization by Xu & Goddard [Xu2004]
 125 → XC_GGA_X_BAYESIAN Bayesian best fit for the enhancement factor [Mortensen2005]
 126 → XC_GGA_X_PBE_JSJR PBE JSJR reparametrization by Pedroza, Silva & Capelle [Pedroza2009]
 130 → XC_GGA_C_PBE Perdew, Burke & Ernzerhof correlation [Perdew1996] [Perdew1997]
 131 → XC_GGA_C_LYP Lee, Yang & Parr [Lee1988] [Miehlich1989]
 132 → XC_GGA_C_P86 Perdew 86 [Perdew1986]
 133 → XC_GGA_C_PBE_SOL Perdew, Burke & Ernzerhof correlation SOL [Perdew2008]
 134 → XC_GGA_C_PW91 Perdew & Wang 91 [Perdew1992]
 135 → XC_GGA_C_AM05 Armiento & Mattsson 05 correlation [Armiento2005] [Mattsson2008]
 136 → XC_GGA_C_XPBE xPBE reparametrization by Xu & Goddard [Xu2004]
 137 → XC_GGA_C_LM Langreth and Mehl correlation [Langreth1981]
 138 → XC_GGA_C_PBE_JRGX JRGX reparametrization by Pedroza, Silva & Capelle [Pedroza2009]
 139 → XC_GGA_X_OPTB88_VDW Becke 88 reoptimized to be used with vdW functional of Dion et al [Klimes2011]
 140 → XC_GGA_X_PBEK1_VDW PBE reparametrization for vdW [Klimes2011]
 141 → XC_GGA_X_OPTPBE_VDW PBE reparametrization for vdW [Klimes2011]
 142 → XC_GGA_X_RGE2 Regularized PBE [Ruzsinszky2009]
 143 → XC_GGA_C_RGE2 Regularized PBE [Ruzsinszky2009]
 144 → XC_GGA_X_RPW86 refitted Perdew & Wang 86 [Murray2009]
 145 → XC_GGA_X_KT1 Keal and Tozer version 1 [Keal2003]
 146 → XC_GGA_XC_KT2 Keal and Tozer version 2 [Keal2003]
Warning
This functional gives NaN on IBM (XG20130608).
 147 → XC_GGA_C_WL Wilson & Levy [Wilson1990]
 148 → XC_GGA_C_WI Wilson & Ivanov [Wilson1998]
 149 → XC_GGA_X_MB88 Modified Becke 88 for proton transfer [Tognetti2009]
 150 → XC_GGA_X_SOGGA Secondorder generalized gradient approximation [Zhao2008]
 151 → XC_GGA_X_SOGGA11 Secondorder generalized gradient approximation 2011 [Peverati2011]
 152 → XC_GGA_C_SOGGA11 Secondorder generalized gradient approximation 2011 [Peverati2011]
 153 → XC_GGA_C_WI0 Wilson & Ivanov initial version [Wilson1998]
 154 → XC_GGA_XC_TH1 Tozer and Handy v. 1 [Tozer1998]
Warning
This functional is not tested. Use at your own risks.
 155 → XC_GGA_XC_TH2 Tozer and Handy v. 2 [Tozer1998a]
 156 → XC_GGA_XC_TH3 Tozer and Handy v. 3 [Handy1998]
 157 → XC_GGA_XC_TH4 Tozer and Handy v. 4 [Handy1998]
 158 → XC_GGA_X_C09X C09x to be used with the VdW of RutgersChalmers [Cooper2010]
 159 → XC_GGA_C_SOGGA11_X To be used with hyb_gga_x_SOGGA11X [Peverati2011a]
 161 → XC_GGA_XC_HCTH_93 HCTH functional fitted to 93 molecules [Hamprecht1998]
 162 → XC_GGA_XC_HCTH_120 HCTH functional fitted to 120 molecules [Boese2000]
 163 → XC_GGA_XC_HCTH_147 HCTH functional fitted to 147 molecules [Boese2000]
 164 → XC_GGA_XC_HCTH_407 HCTH functional fitted to 407 molecules [Boese2001]
 165 → XC_GGA_XC_EDF1 Empirical functionals from Adamson, Gill, and Pople [Adamson1998]
 166 → XC_GGA_XC_XLYP XLYP functional [Xu2004a]
 167 → XC_GGA_XC_B97 Becke 97 [Becke1997]
 168 → XC_GGA_XC_B97_1 Becke 971 [Hamprecht1998] [Becke1997]
 169 → XC_GGA_XC_B97_2 Becke 972 [Becke1997]
 170 → XC_GGA_XC_B97_D Grimme functional to be used with C6 vdW term [Grimme2006]
 171 → XC_GGA_XC_B97_K BoeseMartin for Kinetics [Boese2004]
 172 → XC_GGA_XC_B97_3 Becke 973 [Keal2005]
 173 → XC_GGA_XC_PBE1W Functionals fitted for water [Dahlke2005]
 174 → XC_GGA_XC_MPWLYP1W Functionals fitted for water [Dahlke2005]
 175 → XC_GGA_XC_PBELYP1W Functionals fitted for water [Dahlke2005]
 176 → XC_GGA_XC_SB98_1a SchmiderBecke 98 parameterization 1a [Schmider1998]
 177 → XC_GGA_XC_SB98_1b SchmiderBecke 98 parameterization 1b [Schmider1998]
 178 → XC_GGA_XC_SB98_1c SchmiderBecke 98 parameterization 1c [Schmider1998]
 179 → XC_GGA_XC_SB98_2a SchmiderBecke 98 parameterization 2a [Schmider1998]
 180 → XC_GGA_XC_SB98_2b SchmiderBecke 98 parameterization 2b [Schmider1998]
 181 → XC_GGA_XC_SB98_2c SchmiderBecke 98 parameterization 2c [Schmider1998]
 183 → XC_GGA_X_OL2 Exchange form based on OuYang and Levy v.2 [Fuentealba1995] [OuYang1991]
 184 → XC_GGA_X_APBE mu fixed from the semiclassical neutral atom [Constantin2011]
 186 → XC_GGA_C_APBE mu fixed from the semiclassical neutral atom [Constantin2011]
 191 → XC_GGA_X_HTBS Haas, Tran, Blaha, and Schwarz [Haas2011]
 192 → XC_GGA_X_AIRY Constantin et al based on the Airy gas [Constantin2009]
 193 → XC_GGA_X_LAG Local Airy Gas [Vitos2000]
 194 → XC_GGA_XC_MOHLYP Functional for organometallic chemistry [Schultz2005]
 195 → XC_GGA_XC_MOHLYP2 Functional for barrier heights [Zheng2009]
 196 → XC_GGA_XC_TH_FL Tozer and Handy v. FL [Tozer1997]
 197 → XC_GGA_XC_TH_FC Tozer and Handy v. FC [Tozer1997]
 198 → XC_GGA_XC_TH_FCFO Tozer and Handy v. FCFO [Tozer1997]
 199 → XC_GGA_XC_TH_FCO Tozer and Handy v. FCO [Tozer1997]
 200 → XC_GGA_C_OPTC Optimized correlation functional of Cohen and Handy [Cohen2001]
 (for MetaGGA and Hybrid functionals, with indices in the 200499 range, see the later sections)
 524 → XC_GGA_X_WPBEH shortrange version of the PBE [Heyd2003]
 525 → XC_GGA_X_HJS_PBE HJS screened exchange PBE version [Henderson2008]
 526 → XC_GGA_X_HJS_PBE_SOL HJS screened exchange PBE_SOL version [Henderson2008]
 527 → XC_GGA_X_HJS_B88 HJS screened exchange B88 version [Henderson2008]
Warning
This functional is not tested. Use at your own risks.
 528 → XC_GGA_X_HJS_B97X HJS screened exchange B97x version [Henderson2008]
 529 → XC_GGA_X_ITYH shortrange recipe for exchange GGA functionals [Iikura2001]
Warning
This functional is not tested. Use at your own risks.
MetaGGA functionals (do not forget to add a minus sign, as discussed above). See [Sun2011] for the formulas.
 202 → XC_MGGA_X_TPSS Tao, Perdew, Staroverov & Scuseria [Tao2003] [Perdew2004]
 203 → XC_MGGA_X_M06L Zhao, Truhlar exchange [Zhao2006] [Zhao2007]
 204 → XC_MGGA_X_GVT4 GVT4 (X part of VSXC) from van Voorhis and Scuseria [Voorhis1998]
 205 → XC_MGGA_X_TAU_HCTH tauHCTH from Boese and Handy [Boese2002]
 207 → XC_MGGA_X_BJ06 Becke & Johnson correction to BeckeRoussel 89 [Becke2006]
Warning
This Vxconly mGGA can only be used with a LDA correlation, typically PerdewWang 92 [Perdew1992a], hence ixc =12208 ..
 208 → XC_MGGA_X_TB09 Tranblaha  correction to Becke & Johnson correction to BeckeRoussel 89 [Tran2009]
Warning
This Vxconly mGGA can only be used with a LDA correlation, typically PerdewWang 92 [Perdew1992a].
 209 → XC_MGGA_X_RPP09 Rasanen, Pittalis, and Proetto correction to Becke & Johnson [Rasanen2010]
Warning
This Vxconly mGGA can only be used with a LDA correlation, typically PerdewWang 92 [Perdew1992a].
 232 → XC_MGGA_C_VSXC VSxc from Van Voorhis and Scuseria (correlation part) [Voorhis1998]
Hybrid functionals (do not forget to add a minus sign, as discussed above).
 402 → XC_HYB_GGA_XC_B3LYP The (in)famous B3LYP [Stephens1994]
 406 → XC_HYB_GGA_XC_PBEH PBEH (PBE0) [Adamo1999] [Ernzerhof1999]
 427 → XC_HYB_GGA_XC_HSE03 The 2003 version of the screened hybrid HSE (this case corresponds to hyb_range_fock=\omega^{HF} = 0.15/\sqrt{2} and hyb_range_dft=\omega^{PBE} = 0.15*(2.0)^{1/3} )
 428 → XC_HYB_GGA_XC_HSE06 The 2006 version of the screened hybrid HSE (this case corresponds to hyb_range_fock=hyb_range_dft=\omega^{HF} = \omega^{PBE} = 0.11) [Heyd2003] [Heyd2006] [Krukau2006]
Warning
(The following section is taken from the LibXC sources. In ABINIT, we stick to the LibXC choice.)
Note that there is an enormous mess in the literature concerning the values of omega in HSE. This is due to an error in the original paper that stated that they had used \omega=0.15. This was in fact not true, and the real value used was \omega^{HF} = 0.15 / \sqrt{2} \sim 0.1061 and \omega^{PBE} = 0.15 * (2.0)^{1/3} \sim 0.1890.
In 2006 Krukau et al [Krukau2006] tried to clarify the situation, called HSE03 the above choice of parameters, and called HSE06 to the functional where \omega^{HF}=\omega^{PBE}. By testing several properties for atoms they reached the conclusion that the best value for \omega=0.11. Of course, codes are just as messy as the papers. In Quantum Espresso HSE06 has the value \omega=0.106. VASP, on the other hand, uses for HSE03 the same value \omega^{HF} = \omega^{PBE} = 0.3 (A^{1}) \sim 0.1587, and for HSE06 \omega^{HF} = \omega^{PBE} = 0.2 (A^{1}) \sim 0.1058.
 456 → XC_HYB_GGA_XC_PBE0_13 PBE0⅓ [Cortona2012]
jdtset¶
Mnemonics: index J for DaTaSETs
Characteristics: NO_MULTI
Mentioned in topic(s): topic_multidtset
Variable type: integer
Dimensions: (ndtset)
Default value: [ {‘start’: 1, ‘stop’: ‘ndtset’}; >
Added in version: before_v9
Test list (click to open). Moderately used, [98/1121] in all abinit tests, [8/156] in abinit tutorials
 libxc: t02.abi, t06.abi, t07.abi …
 mpiio: t22.abi …
 paral: t06.abi, t06.abi, t06.abi …
 seq: tsv4_55.abi, tsv4_78.abi, tsv4_80.abi …
 tutoparal: tdmft_1.abi, tdmft_2.abi, tucalc_crpa_1.abi …
 tutoplugs: tw90_5.abi …
 tutorespfn: tpolarization_1.abi, tpolarization_6.abi, tnlo_1.abi …
 v1: t59.abi, t63.abi, t88.abi …
 v2: t26.abi, t75.abi …
 v3: t10.abi, t14.abi, t83.abi …
 v4: t52.abi, t72.abi, t75.abi …
 v5: t03.abi, t21.abi, t55.abi …
 v6: t20.abi, t43.abi, t45.abi …
 v67mbpt: t06.abi, t22.abi, t36.abi …
 v7: t21.abi, t23.abi, t24.abi …
 v8: t01.abi, t01_triqs2_0.abi, t07.abi …
Gives the dataset index of each of the datasets. This index will be used:
 to determine which input variables are specific to each dataset, since the variable names for this dataset will be made from the bare variable name concatenated with this index, and only if such a composite variable name does not exist, the code will consider the bare variable name, or even, the Default;
 to characterize output variable names, if their content differs from dataset to dataset;
 to characterize output files ( root names appended with _DSx where ‘x’ is the dataset index ).
The allowed index values are between 1 and 9999. An input variable name appended with 0 is not allowed. When ndtset == 0, this array is not used, and moreover, no input variable name appended with a digit is allowed. This array might be initialized thanks to the use of the input variable udtset. In this case, jdtset cannot be used.
kpt¶
Mnemonics: K  PoinTs
Mentioned in topic(s): topic_kpoints
Variable type: real
Dimensions: (3,nkpt)
Default value: [0, 0, 0]
Comment: Adequate for one molecule in a supercell
Added in version: before_v9
Test list (click to open). Very frequently used, [1111/1121] in all abinit tests, [156/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t03.abi, t04.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Contains the k points in terms of reciprocal space primitive translations (NOT in cartesian coordinates!). Needed ONLY if kptopt = 0, otherwise deduced from other input variables.
It contains dimensionless numbers in terms of which the cartesian coordinates would be: k_cartesian = k1*G1+k2*G2+k3*G3 where (k1,k2,k3) represent the dimensionless “reduced coordinates” and G1, G2, G3 are the cartesian coordinates of the primitive translation vectors. G1,G2,G3 are related to the choice of direct space primitive translation vectors made in rprim. Note that an overall norm for the k points is supplied by kptnrm. This allows one to avoid supplying many digits for the k points to represent such points as (1,1,1)/3. Note: one of the algorithms used to set up the sphere of G vectors for the basis needs components of kpoints in the range [1,1], so the remapping is easily done by adding or subtracting 1 from each component until it is in the range [1,1]. That is, given the k point normalization kptnrm described below, each component must lie in [ kptnrm, kptnrm ]. Note: a global shift can be provided by %qptn Not read if kptopt/=0.
kptnrm¶
Mnemonics: K  PoinTs NoRMalization
Mentioned in topic(s): topic_kpoints
Variable type: real
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [106/1121] in all abinit tests, [0/156] in abinit tutorials
 builtin: testin_v1.abi, testin_etsf_io.abi …
 etsf_io: t00.abi, t09.abi …
 fast: t03.abi, t05.abi, t06.abi …
 libxc: t03.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 seq: tsv4_90.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t37.abi, t42.abi, t47.abi …
 v3: t01.abi, t52.abi, t54.abi …
 v4: t01.abi, t03.abi, t04.abi …
 v5: t79.abi …
 v7: t60.abi, t61.abi, t62.abi …
Establishes a normalizing denominator for each k point. Needed only if kptopt<=0, otherwise deduced from other input variables. The k point coordinates as fractions of reciprocal lattice translations are therefore kpt(mu,ikpt)/ kptnrm . kptnrm defaults to 1 and can be ignored by the user. It is introduced to avoid the need for many digits in representing numbers such as ⅓. It cannot be smaller than 1.0
kptopt¶
Mnemonics: KPoinTs OPTion
Mentioned in topic(s): topic_kpoints, topic_ElecBandStructure
Variable type: integer
Dimensions: scalar
Default value: 4 if nspden == 4,
1 otherwise.
Added in version: before_v9
Test list (click to open). Very frequently used, [909/1121] in all abinit tests, [96/156] in abinit tutorials
 bigdft: t01.abi, t02.abi, t03.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t02.abi, t09.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_02.abi, tdfpt_03.abi, tdfpt_04.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t03.abi, t05.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Controls the set up of the kpoints list. The aim will be to initialize, by straight reading or by a preprocessing approach based on other input variables, the following input variables, giving the k points, their number, and their weight: kpt, kptnrm, nkpt, and, for iscf/=2, wtk.
Often, the k points will form a lattice in reciprocal space. In this case, one will also aim at initializing input variables that give the reciprocal of this kpoint lattice, as well as its shift with respect to the origin: ngkpt or kptrlatt, as well as on nshiftk and shiftk.
A global additional shift can be provided by %qptn
The use of symmetries (spatial and/or timereversal) is crucial to determine the action of kptopt .
 0 → read directly nkpt, kpt, kptnrm and wtk.

1 → rely on ngkpt or kptrlatt, as well as on nshiftk and shiftk to set up the k points. Take fully into account the symmetry to generate the k points in the Irreducible Brillouin Zone only, with the appropriate weights. (This is the usual mode for GS calculations)

2 → rely on ngkpt or kptrlatt, as well as on nshiftk and shiftk to set up the k points. Take into account only the timereversal symmetry: k points will be generated in half the Brillouin zone, with the appropriate weights. (This is the usual mode when preparing or executing a RF calculation at q=(0 0 0) without noncollinear magnetism)

3 → rely on ngkpt or kptrlatt, as well as on nshiftk and shiftk to set up the k points. Do not take into account any symmetry: k points will be generated in the full Brillouin zone, with the appropriate weights. (This is the usual mode when preparing or executing a RF calculation at nonzero q, or with noncollinear magnetism)

4 → rely on ngkpt or kptrlatt, as well as on nshiftk and shiftk to set up the k points. Take into account all the symmetries EXCEPT the timereversal symmetry to generate the k points in the Irreducible Brillouin Zone, with the appropriate weights. This has to be used when performing calculations with noncollinear magnetism allowed (nspden = 4)

A negative value → rely on kptbounds, and ndivsm (or ndivk) to set up a band structure calculation along different lines (allowed only for iscf == 2). The absolute value of kptopt gives the number of segments of the band structure. Weights are usually irrelevant with this option, and will be left to their default value.
In the case of a grid of k points, the auxiliary variables kptrlen, ngkpt and prtkpt might help you to select the optimal grid.
natom¶
Mnemonics: Number of ATOMs
Mentioned in topic(s): topic_crystal, topic_SmartSymm
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Very frequently used, [1096/1121] in all abinit tests, [144/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives the total number of atoms in the unit cell. Default is 1 but you will obviously want to input this value explicitly. Note that natom refers to all atoms in the unit cell, not only to the irreducible set of atoms in the unit cell (using symmetry operations, this set allows one to recover all atoms). If you want to specify only the irreducible set of atoms, use the symmetriser, see the input variable natrd.
nband¶
Mnemonics: Number of BANDs
Mentioned in topic(s): topic_BandOcc, topic_GW, topic_RPACorrEn, topic_Susceptibility
Variable type: integer
Dimensions: scalar
Default value: None
Comment: the estimated number of occupied bands +1 (TODO provide the mathematical formulation)
Added in version: before_v9
Test list (click to open). Very frequently used, [972/1121] in all abinit tests, [129/156] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t02.abi, t04.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_2.abi, tbase1_3.abi, tbase1_5.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t02.abi, t03.abi, t05.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives number of bands, occupied plus possibly unoccupied, for which wavefunctions are being computed along with eigenvalues. Note: if the parameter occopt (see below) is not set to 2, nband is a scalar integer, but if the parameter occopt is set to 2, then nband must be an array nband (nkpt* nsppol) giving the number of bands explicitly for each k point. This option is provided in order to allow the number of bands treated to vary from k point to k point. For the values of occopt not equal to 0 or 2, nband can be omitted. The number of bands will be set up thanks to the use of the variable fband. The present Default will not be used.
If nspinor is 2, nband must be even for each k point.
In the case of a GW calculation (optdriver = 3 or 4), nband gives the number of bands to be treated to generate the screening (susceptibility and dielectric matrix), as well as the selfenergy. However, to generate the _KSS file (see kssform) the relevant number of bands is given by nbandkss.
nbandhf¶
Mnemonics: Number of BANDs for (Hartree)Fock exact exchange
Mentioned in topic(s): topic_Hybrids
Variable type: integer
Dimensions: scalar
Default value: None
Comment: the estimated number of occupied bands (TODO: provide the mathematical formulation)
Added in version: before_v9
Test list (click to open). Rarely used, [8/1121] in all abinit tests, [0/156] in abinit tutorials
Gives the maximum number of occupied bands with which Fock exact exchange is being computed for the wavefunctions.
ndtset¶
Mnemonics: Number of DaTaSETs
Characteristics: NO_MULTI
Mentioned in topic(s): topic_multidtset
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Very frequently used, [757/1121] in all abinit tests, [88/156] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t11.abi, t20.abi, t21.abi …
 bigdft_paral: t01.abi, t01.abi …
 etsf_io: t02.abi, t21.abi, t30.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t01.abi, t02.abi, t03.abi …
 mpiio: t21.abi, t22.abi, t26.abi …
 paral: t05.abi, t05.abi, t05.abi …
 psml: t02.abi, t03.abi, t04.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdmft_1.abi, tdmft_2.abi, tmbt_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_2.abi, tbase2_1.abi, tbase2_2.abi …
 v1: t25.abi, t39.abi, t43.abi …
 v2: t01.abi, t04.abi, t05.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t01.abi, t02.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t03.abi, t05.abi, t06.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t03.abi, t11.abi, t12.abi …
Gives the number of data sets to be treated. If 0, means that the multidata set treatment is not used, so that the root filenames will not be appended with _DSx, where ‘x’ is the dataset index defined by the input variable jdtset, and also that input names with a dataset index are not allowed. Otherwise, ndtset = 0 is equivalent to ndtset = 1.
ngkpt¶
Mnemonics: Number of Grid points for K PoinTs generation
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: (3)
Default value: [0, 0, 0]
Only relevant if: kptopt >=0,
The use of this variable forbids the use of: kptrlatt
Added in version: before_v9
Test list (click to open). Very frequently used, [668/1121] in all abinit tests, [123/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t07.abi, t14.abi …
 builtin: testin_bigdft.abi …
 etsf_io: t01.abi, t02.abi, t04.abi …
 gpu: t01.abi, t02.abi, t04.abi …
 libxc: t09.abi, t10.abi, t19.abi …
 mpiio: t49.abi, t62.abi, t62.abi …
 paral: t06.abi, t06.abi, t06.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase3_1.abi, tbase3_2.abi, tbase3_3.abi …
 v1: t96.abi …
 v2: t11.abi, t26.abi, t30.abi …
 v3: t01.abi, t13.abi, t14.abi …
 v4: t02.abi, t08.abi, t09.abi …
 v5: t01.abi, t03.abi, t05.abi …
 v6: t04.abi, t06.abi, t08.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t04.abi, t05.abi …
 v8: t02.abi, t04.abi, t05.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Used when kptopt >= 0, if kptrlatt has not been defined (kptrlatt and ngkpt are exclusive of each other). Its three positive components give the number of k points of MonkhorstPack grids (defined with respect to primitive axis in reciprocal space) in each of the three dimensions. ngkpt will be used to generate the corresponding kptrlatt input variable. The use of nshiftk and shiftk, allows one to generate shifted grids, or MonkhorstPack grids defined with respect to conventional unit cells.
When nshiftk = 1, kptrlatt is initialized as a diagonal (3x3) matrix, whose diagonal elements are the three values ngkpt (1:3). When nshiftk is greater than 1, ABINIT will try to generate kptrlatt on the basis of the primitive vectors of the klattice: the number of shifts might be reduced, in which case kptrlatt will not be diagonal anymore.
MonkhorstPack grids are usually the most efficient when their defining integer numbers are even. For a measure of the efficiency, see the input variable kptrlen.
nkpath¶
Mnemonics: Number of Kpoints defining the PATH
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Rarely used, [2/1121] in all abinit tests, [0/156] in abinit tutorials
This variable is used to define the number of highsymmetry kpoints in the kptbounds array when kptopt > 0. Historically, kptbounds is used in conjuction with a negative value of kptopt when performing a NSCF band structure calculation. In this case, the number of kpoints in kptbounds is given by abs(kptopt) + 1. There are, however, other cases in which one has to specify a kpath in the input file in order to activate some kind of post processing tool. Typical examples are the interpolation of the GW corrections at the end of the sigma run or the interpolation of the KS eigenvalues along a path at the end of the SCF run (see also einterp) In a nutshell, nkpath replaces kptopt when we are not performing a NSCF calculation. Note that, unlike kptopt, nkpath represents the total number of points in the kptbounds array.
nkpt¶
Mnemonics: Number of K  Points
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: scalar
Default value: 1 if kptopt == 0,
0 otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [527/1121] in all abinit tests, [41/156] in abinit tutorials
 bigdft: t01.abi, t02.abi, t03.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t04.abi, t09.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 tutoparal: tparal_bandpw_01.abi, tparal_bandpw_02.abi, tparal_bandpw_03.abi …
 tutorespfn: teph4zpr_4.abi, teph4zpr_5.abi, teph4zpr_6.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t07.abi, t11.abi, t23.abi …
 v8: t18.abi, t19.abi, t20.abi …
 v9: t31.abi, t32.abi, t33.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
If nonzero, nkpt gives the number of k points in the k point array kpt. These points are used either to sample the Brillouin zone, or to build a band structure along specified lines.
If nkpt is zero, the code deduces from other input variables (see the list in the description of kptopt) the number of k points, which is possible only when kptopt/=0. If kptopt/=0 and the input value of nkpt /=0, then ABINIT will check that the number of k points generated from the other input variables is exactly the same than nkpt .
If kptopt is positive, nkpt must be coherent with the values of kptrlatt, nshiftk and shiftk. For ground state calculations, one should select the k point in the irreducible Brillouin Zone (obtained by taking into account point symmetries and the timereversal symmetry). For response function calculations, one should select k points in the full Brillouin zone, if the wavevector of the perturbation does not vanish, or in a half of the Brillouin Zone if q=0. The code will automatically decrease the number of k points to the minimal set needed for each particular perturbation.
If kptopt is negative, nkpt will be the sum of the number of points on the different lines of the band structure. For example, if kptopt = 3, one will have three segments; supposing ndivk is 10 12 17, the total number of k points of the circuit will be 10+12+17+1(for the final point)=40.
nkpthf¶
Mnemonics: Number of K  Points for (Hartree) Fock exact exchange
Mentioned in topic(s): topic_Hybrids
Variable type: integer
Dimensions: scalar
Default value: None
Added in version: before_v9
Test list (click to open). Rarely used, [8/1121] in all abinit tests, [0/156] in abinit tutorials
nkpthf gives the number of k points used to sample the full Brillouin zone for the Fock exact exchange contribution. It is obtained from the specification of the wavefunction k point grid (see kptopt), possibly downfolded as specified by fockdownsampling.
nshiftk¶
Mnemonics: Number of SHIFTs for K point grids
Mentioned in topic(s): topic_kpoints
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [559/1121] in all abinit tests, [110/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t07.abi, t14.abi …
 builtin: testin_bigdft.abi, testin_wannier90.abi …
 etsf_io: t02.abi, t04.abi, t30.abi …
 gpu: t02.abi …
 libxc: t09.abi, t19.abi, t41.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t06.abi, t06.abi, t06.abi …
 psml: t03.abi, t04.abi, t07.abi …
 seq: tsv3_03.abi, tsv3_04.abi, tsv4_55.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdmft_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase3_1.abi, tbase3_2.abi, tbase3_3.abi …
 v2: t12.abi, t26.abi, t30.abi …
 v3: t16.abi, t30.abi, t31.abi …
 v4: t02.abi, t08.abi, t09.abi …
 v5: t05.abi, t06.abi, t11.abi …
 v6: t06.abi, t08.abi, t09.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t07.abi, t08.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
This parameter gives the number of shifted grids to be used concurrently to generate the full grid of k points. It can be used with primitive grids defined either from ngkpt or kptrlatt. The maximum allowed value of nshiftk is 8. The values of the shifts are given by shiftk.
The use of nshiftk =1, 2, or 4 is quite common, see the values suggested in the description of shiftk. The other values are either for debugging purposes by experts, or can indicate an error. Such other values are allowed only if chksymbreak=0.
nsppol¶
Mnemonics: Number of SPin POLarization
Mentioned in topic(s): topic_spinpolarisation
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Moderately used, [223/1121] in all abinit tests, [22/156] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t12.abi …
 builtin: testin_v5.abi, testin_libxc.abi …
 etsf_io: t21.abi …
 fast: t17.abi, t19.abi, t20.abi …
 gpu: t03.abi, t05.abi, t05.abi …
 libxc: t00.abi, t02.abi, t05.abi …
 mpiio: t22.abi, t25.abi …
 paral: t05.abi, t05.abi, t05.abi …
 psml: t02.abi, t06.abi, t10.abi …
 seq: tsv2_81.abi, tsv3_03.abi, tsv4_90.abi …
 tutoparal: tpsic_01.abi, tpsic_02.abi, tpsic_03.abi …
 tutorial: tbase1_5.abi, tbase2_1.abi, tbase2_2.abi …
 v1: t08.abi, t09.abi, t21.abi …
 v2: t40.abi, t49.abi, t50.abi …
 v3: t10.abi, t11.abi, t12.abi …
 v4: t08.abi, t59.abi, t79.abi …
 v5: t00.abi, t06.abi, t08.abi …
 v6: t07.abi, t17.abi, t31.abi …
 v67mbpt: t09.abi, t29.abi …
 v7: t05.abi, t06.abi, t07.abi …
 v8: t01.abi, t01_triqs2_0.abi, t18.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t04.abi …
Give the number of INDEPENDENT spin polarisations, for which there are non related wavefunctions. Can take the values 1 or 2.
If nsppol = 1, one has an unpolarized calculation (nspinor = 1, nspden = 1) or an antiferromagnetic system (nspinor = 1, nspden = 2), or a calculation in which spin up and spin down cannot be disentangled (nspinor = 2), that is, either noncollinear magnetism or presence of spin orbit coupling, for which one needs spinor wavefunctions.
If nsppol = 2, one has a spinpolarized (collinear) calculation with separate and different wavefunctions for up and down spin electrons for each band and k point. Compatible only with nspinor == 1, nspden == 2. If nsppol = 2, one usually uses a metallic value for occopt, in order to let ABINIT find the magnetization. On the contrary, if occopt == 1 is used, the user has to impose the magnetization, using spinmagntarget, except for the case of a single isolated Hydrogen atom.
In the present status of development, with nsppol = 1, all values of ixc are allowed, while with nsppol = 2, some values of ixc might not be allowed (e.g. 2, 3, 4, 5, 6, 20, 21, 22 are not allowed).
See also the input variable nspden for the components of the density matrix with respect to the spinpolarization.
nstep¶
Mnemonics: Number of (non)selfconsistent field STEPS
Mentioned in topic(s): topic_SCFControl
Variable type: integer
Dimensions: scalar
Default value: 30
Added in version: before_v9
Test list (click to open). Very frequently used, [1047/1121] in all abinit tests, [135/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_2.abi, teph_tdep_legacy_4.abi, telast_1.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t02.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives the maximum number of cycles (or “iterations”) in a SCF or nonSCF run. Full convergence from random numbers is usually achieved in 1220 SCF iterations. Each can take from minutes to hours. In certain difficult cases, usually related to a small or zero band gap or magnetism, convergence performance may be much worse. When the convergence tolerance tolwfr on the wavefunctions is satisfied, iterations will stop, so for well converged calculations you should set nstep to a value larger than you think will be needed for full convergence, e.g. if using 20 steps usually converges the system, set nstep to 30. For nonselfconsistent runs ( iscf < 0) nstep governs the number of cycles of convergence for the wavefunctions for a fixed density and Hamiltonian.
NOTE that a choice of nstep = 0 is permitted; this will either read wavefunctions from disk (with irdwfk = 1 or irdwfq = 1, or nonzero getwfk or getwfq in the case of multidataset) and compute the density, the total energy and stop, or else (with all of the above vanishing) will initialize randomly the wavefunctions and compute the resulting density and total energy. This is provided for testing purposes. Also NOTE that nstep = 0 with irdwfk = 1 will exactly give the same result as the previous run only if the latter is done with iscf<10 (potential mixing). One can output the density by using prtden. The forces and stress tensor are computed with nstep = 0.
nsym¶
Mnemonics: Number of SYMmetry operations
Mentioned in topic(s): topic_crystal, topic_GW
Variable type: integer
Dimensions: scalar
Default value: 0
Added in version: before_v9
Test list (click to open). Moderately used, [367/1121] in all abinit tests, [19/156] in abinit tutorials
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_bigdft.abi …
 etsf_io: t00.abi, t09.abi, t21.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t04.abi …
 libxc: t01.abi, t13.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 seq: tsv4_80.abi …
 tutoparal: tgswvl_2.abi, tgswvl_2.abi, tgswvl_2.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tpaw1_1.abi, tpaw1_2.abi, tpaw1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t03.abi, t04.abi …
 v5: t01.abi, t06.abi, t09.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v7: t08.abi, t09.abi, t23.abi …
 v8: t05.abi, t67.abi, t68.abi …
 v9: t02.abi, t44.abi, t100.abi …
 vdwxc: t10.abi …
Gives number of space group symmetries to be applied in this problem. Symmetries will be input in array “symrel” and (nonsymmorphic) translations vectors will be input in array “tnons”. If there is no symmetry in the problem then set nsym to 1, because the identity is still a symmetry. In case of a RF calculation, the code is able to use the symmetries of the system to decrease the number of perturbations to be calculated, and to decrease of the number of special k points to be used for the sampling of the Brillouin zone. After the response to the perturbations have been calculated, the symmetries are used to generate as many as possible elements of the 2DTE from those already computed.
Symmetry finder mode (Default mode). If nsym is 0, all the atomic coordinates must be explicitly given (one cannot use the atom manipulator neither the smart symmetrizer): the code will then find automatically the symmetry operations that leave the lattice and each atomic sublattice invariant. It also checks whether the cell is primitive (see chkprim). Note that the tolerance on symmetric atomic positions and lattice is rather stringent: for a symmetry operation to be admitted, the lattice and atomic positions must map on themselves within 1.0e8.
The user is allowed to set up systems with nonprimitive unit cells (i.e. conventional FCC or BCC cells, or supercells without any distortion). In this case, pure translations will be identified as symmetries of the system by the symmetry finder. Then, the combined “pure translation + usual rotation and inversion” symmetry operations can be very numerous. For example, a conventional FCC cell has 192 symmetry operations, instead of the 48 ones of the primitive cell. A maximum limit of 384 symmetry operations is hardcoded. This corresponds to the maximum number of symmetry operations of a 2x2x2 undistorted supercell. Going beyond that number will make the code stop very rapidly. If you want nevertheless, for testing purposes, to treat a larger number of symmetries, change maxnsym.
For GW calculation, the user might want to select only the symmetry operations whose nonsymmorphic translation vector tnons is zero. This can be done with the help of the input variable symmorphi
ntypat¶
Mnemonics: Number of TYPes of AToms
Characteristics: NO_MULTI
Mentioned in topic(s): topic_AtomTypes, topic_PseudosPAW, topic_crystal
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Very frequently used, [1086/1121] in all abinit tests, [144/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v5.abi, testin_bigdft.abi …
 etsf_io: t00.abi, t01.abi, t02.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t01.abi, t02.abi, t03.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t02.abi, t06.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t02.abi, t05.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Gives the number of types of atoms in the unit cell. E.g. for a homopolar system (e.g. pure Si) ntypat is 1.
The code tries to read the same number of pseudopotential files. The first pseudopotential is assigned type number 1, and so on…
There is an exception in the case of alchemical mixing of potentials. In this case, the number of atomic types will differ from the number of pseudopotentials. See mixalch, npsp, ntypalch and %npspalch.
occopt¶
Mnemonics: OCCupation OPTion
Mentioned in topic(s): topic_BandOcc, topic_STM
Variable type: integer
Dimensions: scalar
Default value: 1
Added in version: before_v9
Test list (click to open). Very frequently used, [596/1121] in all abinit tests, [74/156] in abinit tutorials
 atompaw: t04.abi …
 bigdft: t11.abi, t12.abi, t14.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t09.abi, t21.abi …
 fast: t00.abi, t02.abi, t03.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t00.abi, t02.abi, t03.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 psml: t02.abi, t03.abi, t04.abi …
 seq: tsv2_81.abi, tsv2_82.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdmft_1.abi …
 tutoplugs: tw90_5.abi …
 tutorespfn: telast_6.abi, teph_legacy_1.abi, teph4mob_5.abi …
 tutorial: tbase1_5.abi, tbase2_1.abi, tbase2_2.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t06.abi, t09.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t03.abi, t06.abi …
 v7: t01.abi, t02.abi, t04.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t04.abi …
Controls how input parameters nband, occ, and wtk are handled. Possible values are from 0 to 9. For gapped materials (semiconductors, molecules, …), occopt =1 is the favourite for most usages. For metallic situations (also molecules with degenerate levels at Fermi energy), occopt =7 is the favourite for most usages, and one needs moreover to control the input variable tsmear. Use occopt =9 for quasiFermi energy calculations of excited states in gapped materials.

occopt = 0: All k points and spins have the same number of bands. All k points have the same occupancies of bands for a given spin (but these occupancies may differ for spin up and spin down  typical for ferromagnetic insulators). nband is given as a single number, and occ(nband * nsppol) is an array of nband * nsppol elements, read in by the code. The k point weights in array wtk(nkpt) are automatically normalized by the code to add to 1. They cannot differ for differing spins.

occopt = 1: Same as occopt = 0, except that the array occ is automatically generated by the code, to give a semiconductor. An error occurs when filling cannot be done with occupation numbers equal to 2 or 0 in each kpoint (nonspinpolarized case), or with occupation numbers equal to 1 or 0 in each kpoint (spinpolarized case). If nsppol = 2 and occopt == 1 is used, the user has to impose the magnetization, using spinmagntarget, except for the case of a single isolated Hydrogen atom.

occopt = 2: k points may optionally have different numbers of bands and different occupancies. nband(nkpt * nsppol) is given explicitly as an array of nkpt * nsppol elements. occ() is given explicitly for all bands at each k point, and possibly for each spin – the total number of elements is the sum of nband(ikpt) over all k points and spins. The k point weights wtk (nkpt) are NOT automatically normalized under this option.

occopt = 3 to 8 : Metallic occupation of levels, using different occupation schemes (see below). The corresponding thermal broadening, or cold smearing, is defined by the input variable tsmear (see below: the variable xx is the energy in Ha, divided by tsmear). Like for occopt = 1, the variable occ is not read. All k points have the same number of bands, nband is given as a single number, read by the code. The k point weights in array wtk(nkpt) are automatically normalized by the code to add to 1. The combination of a broadening and a physical temperature can be obtained by using both tsmear and tphysel.

occopt = 3: FermiDirac smearing (finitetemperature metal). Smeared delta function: \tilde{\delta}(x)=0.25 (\cosh(x/2.0))^{2}. For usual calculations, at zero temperature, do not use occopt =3, but likely occopt =7. If you want to do a calculation at finite temperature, please also read the information about tphysel.

occopt = 4: “Cold smearing” of N. Marzari (see his thesis work), with a=.5634 (minimization of the bump). Smeared delta function: \tilde{\delta}(x)= (1.5+x(1.5a+x(1.0+ax))) \exp(x^2)/\sqrt{\pi} . Must be used with caution, see the note below.

occopt = 5: “Cold smearing” of N. Marzari (see his thesis work), with a=.8165 (monotonic function in the tail) Same smeared delta function as occopt = 4, with different a. Must be used with caution, see the note below.

occopt = 6: Smearing of Methfessel and Paxton [Methfessel1989] with Hermite polynomial of degree 2, corresponding to “Cold smearing” of N. Marzari with a=0 (so, same smeared delta function as occopt = 4, with different a). Must be used with caution, see the note below.

occopt = 7: Gaussian smearing, corresponding to the 0order Hermite polynomial of Methfessel and Paxton. Smeared delta function: \tilde{\delta}(x)=\exp(x^2)/\sqrt{\pi} . Robust and quite efficient.

occopt = 8: Uniform smearing (the delta function is replaced by a constant function of value one over ]½,½[ (with onehalf value at the boundaries). Used for testing purposes only.

occopt = 9: FermiDirac occupation is enforced with two distinct quasiFermi levels: nqfd holes are forced in bands 1 to ivalence and nqfd electrons are forced in bands with index > ivalence. See details in [Paillard2019]. At present, the number of holes and electrons should be the same. Note that occopt = 9 cannot be used with fixed magnetization calculation.

Note
One can use metallic occupation of levels in the case of a molecule, in order to avoid any problem with degenerate levels. However, it is advised NOT to use occopt = 6 (and to a lesser extent occopt = 4 and 5), since the associated number of electron versus the Fermi energy is NOT guaranteed to be a monotonic function. For true metals, AND a sufficiently dense sampling of the Brillouin zone, this should not happen, but be cautious ! As an indication of this problem, a small variation of input parameters might lead to a jump of total energy, because there might be two or even three possible values of the Fermi energy, and the bisection algorithm finds one or the other.
rprim¶
Mnemonics: Real space PRIMitive translations
Characteristics: EVOLVING
Mentioned in topic(s): topic_UnitCell
Variable type: real
Dimensions: (3,3)
Commentdims: Internally, it is represented as rprim(3,3,nimage)
Default value: [[1, 0, 0], [0, 1, 0], [0, 0, 1]]
Added in version: before_v9
Test list (click to open). Very frequently used, [885/1121] in all abinit tests, [128/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t01.abi, t02.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_v5.abi …
 etsf_io: t00.abi, t02.abi, t04.abi …
 fast: t00.abi, t01.abi, t02.abi …
 gpu: t01.abi, t02.abi, t03.abi …
 libxc: t01.abi, t08.abi, t09.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t02.abi, t03.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdmft_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t01.abi, t06.abi, t07.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t01.abi, t03.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t04.abi, t05.abi, t07.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
Give the three dimensionless primitive translations in real space, to be rescaled by acell and scalecart. The three first numbers are the coordinates of the first vector, the next three numbers are the coordinates of the second, and the last three the coordinates of the third. It is EVOLVING only if ionmov == 2 or 22 and optcell/=0, otherwise it is fixed.
If the Default is used, that is, rprim is the unity matrix, the three dimensionless primitive vectors are three unit vectors in cartesian coordinates. The coordinates (and hence the length) of each vector will be (possibly) multiplied by the corresponding acell value, then (possibly) stretched along the cartesian coordinates by the corresponding scalecart value, to give the dimensional primitive vectors, called %rprimd.
In the general case, the dimensional cartesian coordinates of the crystal primitive translations R1p, R2p and R3p, see %rprimd, are
 R1p(i) = scalecart(i) x rprim (i,1) x acell(1)
 R2p(i) = scalecart(i) x rprim (i,2) x acell(2)
 R3p(i) = scalecart(i) x rprim (i,3) x acell(3)
where i=1,2,3 is the component of the primitive translation (i.e. x, y, and z).
The rprim variable, scaled by scalecart, is thus used to define directions of the primitive vectors, that will be multiplied (so keeping the direction unchanged) by the appropriate length scale acell(1), acell(2), or acell(3), respectively to give the dimensional primitive translations in real space in cartesian coordinates. Presently, it is requested that the mixed product (R1xR2).R3 is positive. If this is not the case, simply exchange a pair of vectors.
To be more specific, keeping the default value of scalecart = 1 to simplify the matter, rprim 1 2 3 4 5 6 7 8 9 corresponds to input of the three primitive translations R1=(1,2,3) (to be multiplied by acell(1)), R2=(4,5,6) (to be multiplied by acell(2)), and R3=(7,8,9) (to be multiplied by acell(3)). Note carefully that the first three numbers input are the first column of rprim , the next three are the second, and the final three are the third. This corresponds with the usual Fortran order for arrays. The matrix whose columns are the reciprocal space primitive translations is the inverse transpose of the matrix whose columns are the direct space primitive translations.
Alternatively to rprim , directions of dimensionless primitive vectors can be specified by using the input variable angdeg. This is especially useful for hexagonal lattices (with 120 or 60 degrees angles). Indeed, in order for symmetries to be recognized, rprim must be symmetric up to tolsym (1.0e5 by default), inducing a specification such as
rprim 0.86602 0.5 0.0
0.86602 0.5 0.0
0.0 0.0 1.0
that can be avoided thanks to angdeg:
angdeg 90 90 120
Note that the following might work as well:
rprim sqrt(0.75) 0.5 0.0
sqrt(0.75) 0.5 0.0
0.0 0.0 1.0
Although the use of scalecart or acell is rather equivalent when the primitive vectors are aligned with the cartesian directions, it is not the case for nonorthogonal primitive vectors. In particular, beginners often make the error of trying to use acell to define primitive vectors in face centered tetragonal lattice, or bodycentered tetragonal lattice, or similarly in face or bodycentered orthorhombic lattices. Let us take the example of a bodycentered tetragonal lattice, that might be defined using the following (“a” and “c” have to be replaced by the appropriate conventional cell vector length):
rprim "a" 0 0
0 "a" 0
"a/2" "a/2" "c/2"
acell 3*1 scalecart 3*1 ! ( These are default values)
The following is a valid, alternative way to define the same primitive vectors:
rprim 1 0 0
0 1 0
1/2 1/2 1/2
scalecart "a" "a" "c"
acell 3*1 ! ( These are default values)
Indeed, the cell has been stretched along the cartesian coordinates, by “a”, “a” and “c” factors.
At variance, the following is WRONG:
rprim 1 0 0
0 1 0
1/2 1/2 1/2
acell "a" "a" "c" ! THIS IS WRONG
scalecart 3*1 ! ( These are default values)
Indeed, the latter would correspond to:
rprim "a" 0 0
0 "a" 0
"c/2" "c/2" "c/2"
acell 3*1 scalecart 3*1 ! ( These are default values)
Namely, the third vector has been rescaled by “c”. It is not at all in the center of the tetragonal cell whose basis vectors are defined by the scaling factor “a”. As another difference between scalecart or acell, note that scalecart is INPUT_ONLY: its content will be immediately applied to rprim, at parsing time, and then scalecart will be set to the default values (3*1). So, in case scalecart is used, the echo of rprim in the output file is not the value contained in the input file, but the value rescaled by scalecart.
rprimd¶
Mnemonics: Real space PRIMitive translations, Dimensional
Characteristics: INTERNAL_ONLY, EVOLVING
Mentioned in topic(s): topic_UnitCell
Variable type: real
Dimensions: (3,3)
Commentdims: Internally, it is represented as rprimd(3,3,nimage).
Default value: None
Added in version: before_v9
This internal variable gives the dimensional real space primitive vectors, computed from acell, scalecart, and rprim.
 R1p(i) = rprimd (i,1) = scalecart(i) x rprim(i,1) x acell(1) for i=1,2,3 (x,y,and z)
 R2p(i) = rprimd (i,2) = scalecart(i) x rprim(i,2) x acell(2) for i=1,2,3
 R3p(i) = rprimd (i,3) = scalecart(i) x rprim(i,3) x acell(3) for i=1,2,3
It is EVOLVING only if ionmov == 2 or 22 and optcell/=0, otherwise it is fixed.
scalecart¶
Mnemonics: SCALE CARTesian coordinates
Characteristics: INPUT_ONLY
Mentioned in topic(s): topic_UnitCell
Variable type: real
Dimensions: (3)
Default value: 3 * 1
Added in version: before_v9
Test list (click to open). Rarely used, [5/1121] in all abinit tests, [0/156] in abinit tutorials
Gives the scaling factors of cartesian coordinates by which dimensionless primitive translations (in “rprim”) are to be multiplied. See the rprim input variable, the acell input variable, and the associated internal %rprimd internal variable. Especially useful for bodycentered and facecentered tetragonal lattices, as well as bodycentered and facecentered orthorhombic lattices, see %rprimd. Note that this input variable is INPUT_ONLY: its content will be immediately applied to rprim, at parsing time, and then scalecart will be set to the default values. So, it will not be echoed.
shiftk¶
Mnemonics: SHIFT for K points
Mentioned in topic(s): topic_kpoints
Variable type: real
Dimensions: (3,nshiftk)
Default value: None if nshiftk>1,
[0.5, 0.5, 0.5] otherwise.
Added in version: before_v9
Test list (click to open). Very frequently used, [624/1121] in all abinit tests, [111/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t00.abi, t07.abi, t14.abi …
 builtin: testin_bigdft.abi, testin_wannier90.abi …
 etsf_io: t02.abi, t04.abi, t30.abi …
 gpu: t02.abi …
 libxc: t09.abi, t19.abi, t41.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t06.abi, t06.abi, t06.abi …
 psml: t01.abi, t03.abi, t04.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdmft_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase3_1.abi, tbase3_2.abi, tbase3_3.abi …
 v1: t96.abi …
 v2: t12.abi, t26.abi, t30.abi …
 v3: t16.abi, t30.abi, t31.abi …
 v4: t02.abi, t08.abi, t09.abi …
 v5: t05.abi, t06.abi, t11.abi …
 v6: t04.abi, t06.abi, t08.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t01.abi, t04.abi, t05.abi …
 v8: t01.abi, t01_triqs2_0.abi, t02.abi …
 v9: t01.abi, t02.abi, t03.abi …
 vdwxc: t10.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
It is used only when kptopt >= 0, and must be defined if nshiftk is larger than 1. shiftk (1:3,1:nshiftk) defines nshiftk shifts of the homogeneous grid of k points based on ngkpt or kptrlatt. The shifts induced by shiftk corresponds to the reduced coordinates in the coordinate system defining the kpoint lattice. For example, if the k point lattice is defined using ngkpt, the point whose reciprocal space reduced coordinates are ( shiftk (1,ii)/ngkpt(1) shiftk (2,ii)/ngkpt(2) shiftk (3,ii)/ngkpt(3) ) belongs to the shifted grid number ii.
The user might rely on ABINIT to suggest suitable and efficient combinations of kptrlatt and shiftk . The procedure to be followed is described with the input variables kptrlen. In what follows, we suggest some interesting values of the shifts, to be used with even values of ngkpt. This list is much less exhaustive than the abovementioned automatic procedure.
1) The default (shifted) MonkhorstPack grids are formed by using nshiftk = 1 and shiftk 0.5 0.5 0.5. This is often the preferred k point sampling, as the shift improves the sampling efficiency with respect to the other simple (nonshifted) possibility nshiftk = 1 and shiftk 0.0 0.0 0.0. There are other interesting possibilities for FCC, BCC and HCP lattices, see later. However, this default can also break symmetry, if the 111 direction is not an axis of rotation. This happens e.g. in tetragonal or hexagonal systems. Abinit will complain about this breaking, and you should adapt shiftk . Easy backup : use nshiftk = 1 and shiftk 0.0 0.0 0.0, to get a nonshifted MonkhorstPack grid, which will be compatible with all symmetries, and is necessary for some features such as kpoint interpolation.
2) When the primitive vectors of the lattice form a FCC lattice, with rprim
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
the (very efficient) usual MonkhorstPack sampling will be generated by using nshiftk = 4 and shiftk
0.5 0.5 0.5
0.5 0.0 0.0
0.0 0.5 0.0
0.0 0.0 0.5
3) When the primitive vectors of the lattice form a BCC lattice, with rprim
0.5 0.5 0.5
0.5 0.5 0.5
0.5 0.5 0.5
the usual MonkhorstPack sampling will be generated by using nshiftk = 2 and shiftk
0.25 0.25 0.25
0.25 0.25 0.25
However, the simple sampling nshiftk = 1 and shiftk 0.5 0.5 0.5 is excellent.
4) For hexagonal lattices with hexagonal axes, e.g. rprim
1.0 0.0 0.0
0.5 sqrt(3)/2 0.0
0.0 0.0 1.0
one can use nshiftk = 1 and shiftk 0.0 0.0 0.5
structure¶
Mnemonics: initialize the crystalline STRUCTURE from …
Mentioned in topic(s): topic_crystal
Variable type: string
Dimensions: scalar
Default value:
Added in version: 9.0.0
Test list (click to open). Moderately used, [294/1121] in all abinit tests, [47/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 builtin: testin_v1.abi …
 etsf_io: t01.abi, t02.abi, t30.abi …
 fast: t03.abi, t05.abi, t06.abi …
 gpu: t02.abi, t03.abi, t04.abi …
 libxc: t19.abi, t41.abi, t42.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t05.abi, t09.abi …
 seq: tsv2_82.abi, tsv3_05.abi, tsv4_78.abi …
 tutoparal: tdfpt_03.abi, tdfpt_04.abi, tdmft_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_4.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase1_1.abi, tbase3_5.abi, tbs_2.abi …
 v1: t00.abi, t38.abi, t59.abi …
 v2: t06.abi, t07.abi, t26.abi …
 v3: t02.abi, t09.abi, t13.abi …
 v4: t03.abi, t04.abi, t05.abi …
 v5: t02.abi, t17.abi, t20.abi …
 v6: t14.abi, t20.abi, t27.abi …
 v67mbpt: t01.abi, t04.abi, t05.abi …
 v7: t07.abi, t10.abi, t22.abi …
 v8: t04.abi, t18.abi, t19.abi …
 v9: t04.abi, t05.abi, t07.abi …
 wannier90: t03.abi …
This variable provides a simplified interface to build the crystalline structure from an external file. The idea is to keep the geometry information separated from the input file so that one can perform multiple calculations in different input files sharing the same structure without having to copy & paste the description of the unit cell inside the input. The single source of truth is now given by an external file that can be easily shared. As a side effect one can easily restart structure relaxations in place by reading the structure from the output file of a previous run.
The structure variable is a string in the format filetype:filepath where:
 filetype specifies the format of the external file
 filepath gives the path to the file relative to the directory where the input file is located.
Variables such as natom, ntypat, typat and znucl are automatically initialized from the external file and need not to be specified in the ABINIT input.
At present (April 25, 2023 ), the allowed values for filetype are:
 abifile → An output file produced by Abinit (only netcdf files are supported for the time being)
 abivars → A txt input file with Abinit variables
 poscar → POSCAR files in VASP5 format (element symbol after the atomic position is required).
Some examples will help clarify.
To read the structure from an external netcdf file produced by Abinit (e.g. out_GSR.nc) use the abifile prefix and the syntax:
structure "abifile:out_GSR.nc"
Other Abinit output files such as the WFK.nc
, the DEN.nc
and the HIST.nc
file
are supported as well e.g.
structure "abifile:out_HIST.nc"
In the case of structural relaxations, these files contain the final geometry (not necessarily relaxed within the given tolerance) hence structure can be used to perform an inplace restart by reading the output of a previous run.
To read the structure from an external file with the structure in Abinit format, use:
structure "abivars:my_text_file"
where my_text_file specifies the lattice in terms [[acell], (%rprimd or angdeg) while the atomic positions are specified with natom and the special variable xred_symbols that can only be used inside this external file.
# MgB2 lattice structure
acell 2*3.086 3.523 Angstrom
rprim 0.866025403784439 0.5 0.0
0.866025403784439 0.5 0.0
0.0 0.0 1.0
natom 3
# Reduced positions followed by element symbol.
xred_symbols
0.0 0.0 0.0 Mg
1/3 2/3 0.5 B
2/3 1/3 0.5 B
To read the structure from an external POSCAR file, use:
structure "poscar:t04_POSCAR"
where t04_POSCAR
is the name of external file.
Critical
Note that the Abinit parser does not support all the possible variants of the POSCAR format. More specifically, we assume a subset of the POSCAR specifications in VASP5 format. This means that the list of element symbols must be present and written on a single line. Only the “direct” and “cartesian” keywords before the atomic positions are allowed (case insensitive).
A typical POSCAR file for hexagonal MgB2 looks like (ignore comments):
Mg1 B2 # Title string (ignored by Abinit)
1.0 # Scaling factor for lattice vectors.
2.672554 1.543000 0.000000 # Lattice vectors in Angstrom (NOTE: ABINIT uses BOHR by default)
2.672554 1.543000 0.000000
0.000000 0.000000 3.523000
Mg B # List of element symbols
1 2 # Number of atoms for each symbol
direct # "direct" for reduced coordinates or "cartesian".
0.000000 0.000000 0.0 Mg # Coordinates followed by the chemical symbol of the atom.
0.333333 0.666667 0.5 B
0.666667 0.333333 0.5 B
A positive scaling factor can be used to rescale the lattice vectors whereas a negative value is interpreted as the volume of the unit cell in Angstrom**3. Obviously, zero is not allowed!
The typat variable is automatically initialized from the list of chemical symbols according to their position in the list. In this example, Mg is of type 1 while B is of type 2.
The ABINIT variables associated to this POSCAR are therefore:
ntypat 2
typat 1 2 2
znucl 12.0 5.0
These are the ONLY QUANTITIES that are initialized from the external POSCAR so please make sure that
your POSCAR resembles the example given above and do not expect ABINIT to understand other entries
such as Selective dynamics
or velocities.
Important
Several POSCAR files available on the internet give atomic positions and lattice vectors with ~6 digits. The ABINIT routines use tighter tolerances to detect the space group thus it may happen that ABINIT does not detect all the symmetry operations with a consequent INCREASE of the number of kpoints in the IBZ and the associated computational cost. This is especially true for hexagonal or rhombohedral lattices. A possible solution is to increase the value of tolsym in the input file to e.g. 1e4 so that ABINIT will automatically refine the atomic positions.
Note the following important remarks:

The structure is initialized by the parser at the very beginning of the calculation hence the external files must exist when Abinit starts to analyze the input file. In a nutshell, structure variables cannot be used to pass the output geometry from one dataset to the next one.

Multidatasets are supported but mind that some variables such as ntypat, typat and znucl are tagged as NO_MULTI. In other words, one can read different files via structure and the multi dataset syntax provided these quantities do not change. ABINIT syntax such as
xred+
are, obviously, not supported. 
The value of typat and znucl given in the input file (if any) is ignored by the parser. The value of natom, ntypat is checked for consistency.

As a rule of thumb, do not try to mix the two approaches: either use structure or the standard (more verbose) approach based on ntypat, typat and znucl to define the unit cell.
Limitations:
 The specification of structures for calculations with images is not supported.
 Alchemical mixing is not supported.
 Reading structure from Fortran file is not yet implemented. It is just a technical problem that will be hopefully solved in the next releases.
In all these cases in which the structure variable is not supported, one has to resort to the standard approach to define the list of atoms and their type.
symrel¶
Mnemonics: SYMmetry in REaL space
Mentioned in topic(s): topic_crystal
Variable type: integer
Dimensions: (3,3,nsym)
Default value: [[1, 0, 0], [0, 1, 0], [0, 0, 1]] if nsym == 1,
None otherwise.
Added in version: before_v9
Test list (click to open). Moderately used, [184/1121] in all abinit tests, [1/156] in abinit tutorials
 bigdft: t00.abi, t01.abi, t02.abi …
 builtin: testin_fast.abi, testin_v1.abi, testin_bigdft.abi …
 fast: t00.abi, t03.abi, t05.abi …
 libxc: t01.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 tutorial: tucalc_lr_2.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t02.abi, t06.abi, t07.abi …
 v4: t03.abi …
 v5: t29.abi, t39.abi, t48.abi …
 v6: t01.abi, t02.abi, t12.abi …
 v7: t32.abi, t60.abi, t61.abi …
 v8: t95.abi …
 v9: t17.abi, t44.abi …
Gives “nsym” 3x3 matrices expressing space group symmetries in terms of their action on the direct (or real) space primitive translations. It turns out that these can always be expressed as integers. Always give the identity matrix even if no other symmetries hold, e.g. symrel 1 0 0 0 1 0 0 0 1.
Also note that for this array, as for all others, the array elements are filled in a columnwise order as is usual for Fortran. Explicitly, symrel 1 0 0 1 1 0 0 0 1 for symmetry operation isym is stored internally as symrel(1,1,isym)=1, symrel(1,2)=1, … The atom located at xred(1:3) is send to location xred_sym(jj)=symrel(jj,1,isym)*xred(1)+symrel(jj,2,isym)*xred(2)+symrel(jj,3,isym)*xred(3)+tnons(jj).
The relation between the above symmetry matrices symrel , expressed in the basis of primitive translations, and the same symmetry matrices expressed in cartesian coordinates, is as follows. Denote the matrix whose columns are the primitive translations as R, and denote the cartesian symmetry matrix as S. Then symrel = R(inverse) * S * R where matrix multiplication is implied. When the symmetry finder is used (see nsym), symrel will be computed automatically. Also see the accompanying input variables tnons and symafm, for the full definition of the symmetry operations. Such variables are used to infer spgroup, spgroupma and genafm if they are not userdefined.
tnons¶
Mnemonics: Translation NONSymmorphic vectors
Mentioned in topic(s): topic_crystal
Variable type: real
Dimensions: (3,nsym)
Default value: None
Added in version: before_v9
Test list (click to open). Moderately used, [213/1121] in all abinit tests, [4/156] in abinit tutorials
 builtin: testin_fast.abi, testin_v1.abi …
 fast: t00.abi, t03.abi, t05.abi …
 libxc: t01.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv4_55.abi …
 tutoparal: tdfpt_03.abi, tdfpt_04.abi …
 tutorial: tbase4_7.abi, tbase4_8.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t06.abi, t07.abi, t08.abi …
 v4: t01.abi, t03.abi, t46.abi …
 v5: t29.abi, t48.abi, t49.abi …
 v6: t01.abi, t02.abi, t12.abi …
 v7: t09.abi, t60.abi, t61.abi …
 v8: t12.abi, t24.abi, t25.abi …
 v9: t01.abi, t02.abi, t03.abi …
Gives the (nonsymmorphic) translation vectors associated with the symmetries expressed in “symrel”. These may all be 0, or may be fractional (nonprimitive) translations expressed relative to the real space primitive translations (so, using the “reduced” system of coordinates, see “xred”). If all elements of the space group leave 0 0 0 invariant, then these are all 0. When the symmetry finder is used (see nsym), tnons is computed automatically.
For the groundstate and DFPT drivers of ABINIT, the value of tnons is unrestricted. However, for GW and BSE, the symmetry operations must leave the FFT grid invariant. Preparatory (Groundstate) runs must also use the same atomic geometry, hence the same tnons. As ABINIT cannot guess whether the user has in mind to do a GW or BSE run after the GS run, a conservative approach is implemented, requiring such match of symmetry operations and FFT grid also in the GS case. See more details in the section describing the input variable chksymtnons.
See also symafm for the complete description of the symmetry operation.
toldfe¶
Mnemonics: TOLerance on the DiFference of total Energy
Characteristics: ENERGY
Mentioned in topic(s): topic_SCFControl
Variable type: real
Dimensions: scalar
Default value: 0.0
Comment: The default value implies that this stopping condition is ignored. For the SCF case, one and only one of the input tolerance criteria tolwfr, toldff, tolrff, toldfe or tolvrs must differ from zero.
The use of this variable forbids the use of: tolwfr or toldff or tolrff or tolvrs
Added in version: before_v9
Test list (click to open). Moderately used, [320/1121] in all abinit tests, [41/156] in abinit tutorials
 bigdft: t00.abi, t01.abi, t02.abi …
 builtin: testin_bigdft.abi, testin_etsf_io.abi …
 etsf_io: t00.abi, t01.abi, t04.abi …
 gpu: t04.abi …
 libxc: t03.abi, t51.abi, t52.abi …
 mpiio: t21.abi, t22.abi, t25.abi …
 paral: t06.abi, t06.abi, t06.abi …
 psml: t02.abi, t03.abi, t04.abi …
 seq: tsv3_03.abi, tsv3_04.abi, tsv4_55.abi …
 tutoparal: tparal_bandpw_01.abi, tparal_bandpw_02.abi, tparal_bandpw_03.abi …
 tutorespfn: tpolarization_1.abi, tpolarization_4.abi, tpolarization_6.abi …
 tutorial: tbase1_1.abi, tbase1_2.abi, tbase1_3.abi …
 v1: t08.abi, t09.abi, t40.abi …
 v2: t11.abi, t36.abi, t43.abi …
 v3: t01.abi, t02.abi, t14.abi …
 v4: t04.abi, t05.abi, t06.abi …
 v5: t01.abi, t07.abi, t08.abi …
 v6: t11.abi, t16.abi, t17.abi …
 v67mbpt: t11.abi, t15.abi, t31.abi …
 v7: t05.abi, t06.abi, t09.abi …
 v8: t02.abi, t03.abi, t12.abi …
 v9: t29.abi, t33.abi, t34.abi …
 vdwxc: t10.abi …
 wannier90: t03.abi …
Sets a tolerance for absolute differences of total energy that, reached TWICE successively, will cause one SCF cycle to stop (and ions to be moved). Can be specified in Ha (the default), Ry, eV or Kelvin, since toldfe has the ENERGY characteristics (1 Ha = 27.2113845 eV). If set to zero, this stopping condition is ignored. Effective only when SCF cycles are done (iscf>0). Because of machine precision, it is not worth to try to obtain differences in energy that are smaller than about 1.0d12 of the total energy. To get accurate stresses may be quite demanding. When the geometry is optimized (relaxation of atomic positions or primitive vectors), the use of toldfe is to be avoided. The use of toldff or tolrff is by far preferable, in order to have a handle on the geometry characteristics. When all forces vanish by symmetry (e.g. optimization of the lattice parameters of a highsymmetry crystal), then place toldfe to 1.0d12, or use (better) tolvrs. Since toldfe , toldff, tolrff, tolvrs and tolwfr are aimed at the same goal (causing the SCF cycle to stop), they are seen as a unique input variable at reading. Hence, it is forbidden that two of these input variables have nonzero values for the same dataset, or generically (for all datasets). However, a nonzero value for one such variable for one dataset will have precedence on the nonzero value for another input variable defined generically.
toldff¶
Mnemonics: TOLerance on the DiFference of Forces
Mentioned in topic(s): topic_SCFControl, topic_ForcesStresses
Variable type: real
Dimensions: scalar
Default value: 0.0
Comment: The default value implies that this stopping condition is ignored. For the SCF case, one and only one of the input tolerance criteria tolwfr, toldff, tolrff, toldfe or tolvrs must differ from zero.
The use of this variable forbids the use of: tolwfr or toldfe or tolrff or tolvrs
Added in version: before_v9
Test list (click to open). Moderately used, [102/1121] in all abinit tests, [15/156] in abinit tutorials
 builtin: testin_fast.abi …
 etsf_io: t21.abi …
 fast: t00.abi …
 libxc: t10.abi …
 mpiio: t21.abi, t22.abi, t24.abi …
 paral: t03.abi, t03.abi, t03.abi …
 seq: tsv6_121.abi, tsv6_122.abi, tsv6_123.abi …
 tutoparal: timages_01.abi, timages_02.abi, timages_03.abi …
 tutorial: tbase1_3.abi, tbase2_1.abi, tbase2_2.abi …
 v1: t41.abi, t44.abi, t45.abi …
 v2: t40.abi, t41.abi, t44.abi …
 v3: t40.abi, t43.abi, t80.abi …
 v4: t42.abi …
 v5: t02.abi, t03.abi, t07.abi …
 v6: t06.abi, t15.abi, t21.abi …
 v7: t08.abi …
 v8: t05.abi …
 v9: t22.abi, t75.abi, t104.abi …
Sets a tolerance for differences of forces (in hartree/Bohr) that, reached TWICE successively, will cause one SCF cycle to stop (and ions to be moved). If set to zero, this stopping condition is ignored. Effective only when SCF cycles are done (iscf>0). This tolerance applies to any particular cartesian component of any atom, INCLUDING fixed ones. This is to be used when trying to equilibrate a structure to its lowest energy configuration (select ionmov), or in case of molecular dynamics (ionmov = 1) A value ten times smaller than tolmxf is suggested (for example 5.0d6 hartree/Bohr). This stopping criterion is not allowed for RF calculations. Since toldfe, toldff , tolrff, tolvrs and tolwfr are aimed at the same goal (causing the SCF cycle to stop), they are seen as a unique input variable at reading. Hence, it is forbidden that two of these input variables have nonzero values for the same dataset, or generically (for all datasets). However, a nonzero value for one such variable for one dataset will have precedence on the nonzero value for another input variable defined generically.
tolrff¶
Mnemonics: TOLerance on the Relative diFference of Forces
Mentioned in topic(s): topic_SCFControl, topic_ForcesStresses
Variable type: real
Dimensions: scalar
Default value: 0.0
Comment: The default value implies that this stopping condition is ignored. For the SCF case, one and only one of the input tolerance criteria tolwfr, toldff, tolrff, toldfe or tolvrs must differ from zero.
The use of this variable forbids the use of: tolwfr or toldfe or toldff or tolvrs‘
Added in version: before_v9
Test list (click to open). Moderately used, [14/1121] in all abinit tests, [3/156] in abinit tutorials
 bigdft: t22.abi
 libxc: t67.abi, t68.abi, t69.abi, t70.abi, t71.abi
 tutoparal: tpsic_02.abi, tpsic_03.abi
 tutorial: tbase1_3.abi
 v5: t41.abi
 v8: t17.abi, t30.abi, t31.abi
 v9: t18.abi
Sets a tolerance for the ratio of differences of forces (in hartree/Bohr) to maximum force, that, reached TWICE successively, will cause one SCF cycle to stop (and ions to be moved): diffor < tolrff * maxfor. If set to zero, this stopping condition is ignored. Effective only when SCF cycles are done (iscf>0). This tolerance applies to any particular cartesian component of any atom, INCLUDING fixed ones. This is to be used when trying to equilibrate a structure to its lowest energy configuration (select ionmov), or in case of molecular dynamics (ionmov = 1) A value of 0.02 is suggested. This stopping criterion is not allowed for RF calculations. Since toldfe, toldff, tolrff , tolvrs and tolwfr are aimed at the same goal (causing the SCF cycle to stop), they are seen as a unique input variable at reading. Hence, it is forbidden that two of these input variables have nonzero values for the same dataset, or generically (for all datasets). However, a nonzero value for one such variable for one dataset will have precedence on the nonzero value for another input variable defined generically.
tolvrs¶
Mnemonics: TOLerance on the potential V(r) ReSidual
Mentioned in topic(s): topic_SCFControl
Variable type: real
Dimensions: scalar
Default value: 0.0
Comment: The default value implies that this stopping condition is ignored. For the SCF case, one and only one of the input tolerance criteria tolwfr, toldff, tolrff, toldfe or tolvrs must differ from zero.
The use of this variable forbids the use of: tolwfr or toldfe or toldff or tolrff‘
Added in version: before_v9
Test list (click to open). Moderately used, [377/1121] in all abinit tests, [67/156] in abinit tutorials
 atompaw: t02.abi, t04.abi …
 bigdft: t23.abi, t31.abi, t32.abi …
 bigdft_paral: t01.abi, t01.abi, t02.abi …
 gpu: t02.abi, t03.abi, t05.abi …
 libxc: t13.abi, t41.abi, t42.abi …
 mpiio: t28.abi, t62.abi, t62.abi …
 paral: t28.abi, t31.abi, t32.abi …
 tutoparal: tdfpt_01.abi, tdfpt_02.abi, tdfpt_03.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbasepar_1.abi, tbasepar_2.abi, tbs_1.abi …
 v1: t08.abi, t39.abi, t53.abi …
 v2: t11.abi, t12.abi, t26.abi …
 v3: t02.abi, t07.abi, t08.abi …
 v4: t09.abi, t17.abi, t20.abi …
 v5: t05.abi, t06.abi, t07.abi …
 v6: t07.abi, t14.abi, t35.abi …
 v67mbpt: t04.abi, t06.abi, t07.abi …
 v7: t01.abi, t02.abi, t03.abi …
 v8: t01.abi, t01_triqs2_0.abi, t04.abi …
 v9: t01.abi, t02.abi, t03.abi …
 wannier90: t11.abi, t12.abi, t13.abi …
Sets a tolerance for potential residual that, when reached, will cause one SCF cycle to stop (and ions to be moved). If set to zero, this stopping condition is ignored. Effective only when SCF cycles are done (iscf>0). To get accurate stresses may be quite demanding. For simple materials with internal positions determined by symmetries, a value of tolvrs = 10^12 empirically leads to a very approximate 10^6 atomic unit accuracy for the optimized lattice parameter.
Additional explanation: the residual of the potential is the difference between the input potential and the output potential, when the latter is obtained from the density determined from the eigenfunctions of the input potential. When the selfconsistency loop is achieved, both input and output potentials must be equal, and the residual of the potential must be zero. The tolerance on the potential residual is imposed by first subtracting the mean of the residual of the potential (or the trace of the potential matrix, if the system is spinpolarized), then summing the square of this function over all FFT grid points. The result should be lower than tolvrs . Since toldfe, toldff, tolrff, tolvrs and tolwfr are aimed at the same goal (causing the SCF cycle to stop), they are seen as a unique input variable at reading. Hence, it is forbidden that two of these input variables have nonzero values for the same dataset, or generically (for all datasets). However, a nonzero value for one such variable for one dataset will have precedence on the nonzero value for another input variable defined generically.
tolwfr¶
Mnemonics: TOLerance on WaveFunction squared Residual
Mentioned in topic(s): topic_SCFControl
Variable type: real
Dimensions: scalar
Default value: 0.0
Comment: The default value implies that this stopping condition is ignored. For the SCF case, one and only one of the input tolerance criteria tolwfr, toldff, tolrff, toldfe or tolvrs must differ from zero.
The use of this variable forbids the use of: toldfe or toldff or tolrff or tolvrs
Added in version: before_v9
Test list (click to open). Moderately used, [527/1121] in all abinit tests, [55/156] in abinit tutorials
 bigdft: t06.abi, t07.abi, t10.abi …
 builtin: testin_v1.abi, testin_v5.abi, testin_libxc.abi …
 etsf_io: t02.abi, t22.abi, t30.abi …
 fast: t01.abi, t02.abi, t03.abi …
 gpu: t01.abi …
 libxc: t00.abi, t01.abi, t02.abi …
 mpiio: t01.abi, t01.abi, t01.abi …
 paral: t01.abi, t01.abi, t01.abi …
 psml: t01.abi, t05.abi, t09.abi …
 seq: tsv2_81.abi, tsv2_82.abi, tsv3_03.abi …
 tutoparal: tdfpt_04.abi, tdmft_1.abi, tgswvl_1.abi …
 tutoplugs: tw90_1.abi, tw90_2.abi, tw90_3.abi …
 tutorespfn: teph_tdep_legacy_1.abi, teph_tdep_legacy_2.abi, teph_tdep_legacy_3.abi …
 tutorial: tbase3_5.abi, tbs_1.abi, tgw1_1.abi …
 v1: t00.abi, t01.abi, t02.abi …
 v2: t01.abi, t02.abi, t03.abi …
 v3: t06.abi, t09.abi, t13.abi …
 v4: t01.abi, t02.abi, t03.abi …
 v5: t00.abi, t05.abi, t14.abi …
 v6: t01.abi, t02.abi, t03.abi …
 v67mbpt: t01.abi, t02.abi, t03.abi …
 v7: t07.abi, t12.abi, t13.abi …
 v8: t04.abi, t07.abi, t24.abi …
 v9: t05.abi, t07.abi, t09.abi …
 wannier90: t00.abi, t01.abi, t02.abi …
The signification of this tolerance depends on the basis set. In plane waves, it gives a convergence tolerance for the largest squared “residual” (defined below) for any given band. The squared residual is: < nk (HE)^2 nk>, E = < nkHnk >
which clearly is nonnegative and goes to 0 as the iterations converge to an eigenstate. With the squared residual expressed in Hartrees^2 (Hartrees squared), the largest squared residual (called residm in the code) encountered over all bands and kpoints must be less than tolwfr for iterations to halt due to successful convergence. Note that if iscf > 0, this criterion should be replaced by those based on toldfe (preferred for ionmov == 0), toldff tolrff (preferred for ionmov /= 0), or tolvrs (preferred for theoretical reasons!). When tolwfr is 0.0, this criterion is ignored, and a finite value of toldfe, toldff or tolvrs must be specified. This also imposes a restriction on taking an ion step; ion steps are not permitted unless the largest squared residual is less than tolwfr, ensuring accurate forces. To get accurate stresses may be quite demanding. Note that the preparatory GS calculations before a RF calculations must be highly converged. Typical values for these preparatory runs are tolwfr between 1.0d16 and 1.0d22.
Note that tolwfr is often used in the test cases, but this is tolwfr purely for historical reasons: except when iscf < 0, other criteria should be used.
In the wavelet case (see usewvl = 1), this criterion is the favoured one. It is based on the norm 2 of the gradient of the wavefunctions. Typical values range from 5*10^4 to 5*10^5.
Since toldfe, toldff, tolrff, tolvrs and tolwfr are aimed at the same goal (causing the SCF cycle to stop), they are seen as a unique input variable at reading. Hence, it is forbidden that two of these input variables have nonzero values for the same dataset, or generically (for all datasets). However, a nonzero value for one such variable for one dataset will have precedence on the nonzero value for another input variable defined generically.
typat¶
Mnemonics: TYPe of AToms
Mentioned in topic(s): topic_crystal, topic_AtomTypes
Variable type: integer
Dimensions: [3, ‘natrd’] if natrd<natom,
[3, ‘natom’] otherwise.
Default value: 1 if natom == 1, None otherwise.
Added in version: before_v9